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The ionization conjecture 
in Hartree-Fock theory 

By Jan Philip Solovej* 
Abstract 

We prove the ionization conjecture within the Hartree-Fock theory of 
atoms. More precisely, we prove that, if the nuclear charge is allowed to tend 
to infinity, the maximal negative ionization charge and the ionization energy of 
atoms nevertheless remain bounded. Moreover, we show that in Hartree-Fock 
theory the radius of an atom (properly defined) is bounded independently of 
its nuclear charge. 
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1. Introduction and main results 



One of the great triumphs of quantum mechanics is that it explains the 
order in the periodic table qualitatively as well as quantitatively. In elementary 
chemistry it is discussed how quantum mechanics implies the shell structure 
of atoms which gives a qualitative understanding of the periodic table. In 
computational quantum chemistry it is found that quantum mechanics gives 
excellent agreement with the quantitative aspects of the periodic table. It is 
a very striking fact, however, that the periodic table is much more "periodic" 
than can be explained by the simple shell structure picture. As an example it 
can be mentioned that e.g., the radii of different atoms belonging to the same 
group in the periodic table do not vary very much, although the number of 
electrons in the atoms can vary by a factor of 10. Another related example is 
the fact that the maximal negative ionization (the number of extra electrons 
that a neutral atom can bind) remains small (possibly no bigger than 2) even 
for atoms with large atomic number (nuclear charge). These experimental facts 
can to some extent be understood numerically, but there is no good qualitative 
explanation for them. 

In the mathematical physics literature the problem has been formulated 
as follows (see e.g., Problems 10C and 10D in [22] or Problems 9 and 10 in 
[23]). Imagine that we consider 'the infinitely large periodic table', i.e., atoms 
with arbitrarily large nuclear charge Z; is it then still true that the radius and 
maximal negative ionization remain bounded? This question often referred to 
as the ionization conjecture is the subject of this paper. 

To be completely honest neither the qualitative nor the quantitative expla- 
nations of the periodic table use the full quantum mechanical description. On 
one hand the simple qualitative shell structure picture ignores the interactions 
between the electrons in the atoms. On the other hand even in computational 
quantum chemistry one most often uses approximations to the full many body 
quantum mechanical description. There are in fact a hierarchy of models for the 
structure of atoms. The one which is usually considered most complete is the 
Schrodinger many-particle model. There are, however, even more complicated 
models, which take relativistic and/or quantum field theoretic corrections into 
account. 

A description which is somewhat simpler than the Schrodinger model is the 
Hartree-Fock (HF) model. Because of its greater simplicity it has been more 
widely used in computational quantum chemistry than the full Schrodinger 
model. Although, chemists over the years have developed numerous gener- 
alizations of the Hartree-Fock model, it is still remarkable how tremendously 
successful the original (HF) model has been in describing the structure of atoms 
and molecules. 
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A model which is again much simpler than the Hartree-Fock model is the 
Thomas- Fermi (TF) model. In this model the problem of finding the structure 
of an atom is essentially reduced to solving an ODE. The TF model has some 
features, which are qualitatively wrong. Most notably it predicts that atoms 
do not bind to form molecules (Teller's no binding theorem; see [17]). 

In this work we shall show that the TF model is, indeed, a much better 
approximation to the more complicated HF model than generally believed. In 
fact, we shall show that it is only the outermost region of the atom which is 
not well described by the TF model. 

As a simple corollary of this improved TF approximation we shall prove 
the ionization conjecture within HF theory. The corresponding results for the 
full Schrodinger theory are still open and only much simpler results are known 
(see e.g., [5], [15], [20], [21], [24]). In [3] the ionization conjecture was solved in 
the Thomas-Fermi- von Weizsacker generalization of the Thomas- Fermi model. 
In [25] the ionization conjecture was solved in a simplified Hartree-Fock mean 
field model by a method very similar to the one presented here. In the simplified 
model the atoms are entirely spherically symmetric. In the full HF model, 
however, the atoms need not be spherically symmetric. This lack of spherical 
symmetry in the HF model is one of the main reasons for many of the difficulties 
that have to be overcome in the present paper, although this may not always 
be apparent from the presentation. 

We shall now describe more precisely the results of this paper. In common 
for all the atomic models is that, given the number of electrons N and the 
nuclear charge Z, they describe how to find the electronic ground state density 
p £ L 1 (IR 3 ), with j p = N. Or more precisely how to find one ground state 
density, since it may not be unique. In the TF model the ground state is 
described only by the density, whereas in the Schrodinger and HF models the 
density is derived from more detailed descriptions of the ground state. For all 
models we shall use the following definitions. We distinguish quantities in the 
different models by adding superscripts TF, HF. (In this work we shall not be 
concerned with the Schrodinger model at all.) Throughout the paper we use 
units in which % = m = e = 1, i.e., atomic units. 

We shall discuss Hartree-Fock theory in greater detail in Section 3 and 
Thomas-Fermi theory in greater detail in Section 4. For a complete discussion 
of TF theory we refer the reader to the original paper by Lieb and Simon [17] 
or the review by Lieb [10]. In this introduction we shall only make the most 
basic definitions and enough remarks in order to state some of the main results 
of the paper. 

Definition 1.1. (Mean field potentials). Let p HF and p TF be the densities 
of atomic ground states in the HF and TF models respectively. We define the 
corresponding mean field potentials 
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(1) (p UF (x) := Z\x\- 1 -p BF *\x\- 1 = Z\x\- 1 - Jp^Wlx-yl^dy 

(2) ip TF {x) := Z\x\- 1 -p TF *\x\- 1 = Z\x\- 1 - Jp^Wlx-yl^dy 
and for all R > the screened nuclear potentials at radius R 

(3) <S>™(x) := Z\x\~ x -\ p^{y)\x-y\^dy 

J\y\<R 

(4) ;= Z\x\^-\ p^{y)\x-y\- x dy. 

J\y\<R 

This is the potential from the nuclear charge Z screened by the electrons in the 
region {x : \x\ < R}. The screened nuclear potential will be very important in 
the technical proofs in SectionslO-13. 

Definition 1.2. (Radius). Let again p HF and p TF be the densities of atomic 
ground states in the HF and TF models respectively. We define the radius 
Rz,n{v) to the v last electrons by 

/ P TF (x) dx = u, I P UF (x) dx = v. 

J\x\>R^ N (u) J\x\>R^ N (u) 

The functions (p TF and p TF are the unique solutions to the set of equations 

(5) A<p TF (x) = 4ttP tf {x)-4ttZ5(x) 

(6) P TF (x) = 2 3 / 2 (3vr 2 )^ 1 [^ TF (x)-M TF ]5 /2 

(7) J P TF = N. 

Here p TF is a nonnegative parameter called the chemical potential, which is 
also uniquely determined from the equations. We have used the notation 
[t]+ = max{t, 0} for all f gl. The equations (5-7) only have solutions when 
N < Z. For N > Z we shall let (p TF and p TF refer to the solutions for N = Z, 
the neutral case. Instead of fixing N and determining p TF (the 'canonical' pic- 
ture) one could fix p TF and determine N (the 'grand canonical' picture). The 
equation (5) is essentially equivalent to (2) and expresses the fact that (p TF is 
the mean field potential generated by the positive charge Z and the negative 
charge distribution —p TF . The equations (6-7) state that p TF is given by the 
semiclassical expression for the density of an electron gas of N electrons in the 
exterior potential ip TF . For a discussion of semiclassics we refer the reader to 
Section 8. 

Remark 1.3. The total energy of the atom in Thomas-Fermi theory is 

(8) ^(3vr 2 ) 2 / 3 / p TF (xf/ 3 dx-Z f p TP {x)\x\- 1 dx 



+\JJp TV {x)\x-y\- l p FF {y)dxdy > -e Z 7 / 3 
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where eo is the total binding energy of a neutral TF atom of unit nuclear 
charge. Numerically [10], 

(9) e = 2(3vr 2 )- 2 / 3 • 3.67874 = 0.7687. 

For a neutral atom, where N = Z, the above inequality is an equality. The 
inequality states that in Thomas-Fermi theory the energy is smallest for a 
neutral atom. 

We can now state two of the main results in this paper. 

Theorem 1.4 (Potential estimate). For all Z > 1 and all integers N 
with N > Z for which there exist Hartree-Fock ground states with J p HF = N 
we have 

(10) |v? HF (x) - v? TF (x)| < ^|x|- 4+£ ° + A u 
where A,p,Ai,£o > are universal constants. 

This theorem is proved in Section 13 on page 535. The significance of the 
power |x| -4 is that for N > Z we have lim,£_ >00 (p TF (x) = 3 4 2~ 3 7r 2 |x|~ 4 . The 
existence of this limit known as the Sommerfeld asymptotic law [27] follows 
from Theorem 2.10 in [10], but we shall also prove it in Theorems 5.2 and 5.4 
below. 

Note that the bound in Theorem 1.4 is uniform in N and Z. 

The second main theorem is the universal bound on the atomic radius 
mentioned in the beginning of the introduction. In fact, not only do we prove 
uniform bounds but we also establish a certain exact asymptotic formula for 
the radius of an "infinite atom" . 

Theorem 1.5. Both liminf R z F z (v) and limsup-R§ F z (^) are bounded 
and have the asymptotic behavior 

2-1/334/3^2/3^-1/3 + 0^-1/3) 

as v — ► 00. 

The proof of this theorem can be found in Section 13 on page 535. The 
universal bound on the maximal ionization is given in Theorem 3.6. The proof 
is given in Section 13 on page 534. A universal bound on the ionization energy 
(the energy it takes to remove one electron) is formulated in Theorem 3.8. 
The proof is given in Section 13 on page 537. Theorems 3.6 and 3.8 are 
as important as Theorems 1.4 and 1.5. We have deferred the statements of 
Theorems 3.6 and 3.8 in order not to have to make too many definitions here 
in the introduction. 

One of the main ideas in the paper is to use the strong universal behav- 
ior of the TF theory reflected in the Sommerfeld asymptotics. If we com- 
bine (5) and (6) we see that for ^ TF = the potential satisfies the equation 
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A(p TF (x) = 2 7 / 2 (37r)~ 1 [(y9 TF ]^ 2 (x) for x / 0. It turns out that the singularity 
at x = of any solution to this equation is either of weak type ~ for 
some constant Z or of strong type ~ 3 4 2~ 3 7r 2 |x| -4 (see [30] for a discussion 
of singularities for differential equations of similar type). The surprising fact, 
contained in Theorem 1.4, is that the same type of universal behavior holds 
also for the much more complicated HF potential. We prove this by comparing 
with appropriately modified TF systems on different scales, using the fact that 
the modifications do not affect the universal behavior. A direct comparison 
works only in a short range of scales. This is however enough to use an iter- 
ative renormalization argument to bootstrap the comparison to essentially all 
scales. 

The paper is organized as follows. In Section 2 we fix our notational 
conventions and give some basic prerequisites. In Section 3 we discuss Hartree- 
Fock theory. In Sections 4 and 5 we discuss Thomas-Fermi theory. In particular 
we show that the TF model, indeed, has the universal behavior for large Z that 
we want to establish for the HF model. In the TF model the universality can 
be expressed very precisely through the Sommerfeld asymptotics. 

In Section 6 we begin the more technical work. We show in this section 
that the HF atom in the region {x : \x\ > R} is determined to a good approx- 
imation, in terms of energy, from knowledge of the screened nuclear potential 
<3? R . It is this crucial step in the whole argument that I do not know how to 
generalize to the Schrodinger model or even to the case of molecules in HF 
theory. 

For the outermost region of the atom one cannot use the energy to control 
the density. In fact, changing the density of the atom far from the nucleus will 
not affect the energy very much. Far away from the nucleus one must use the 
exact energy minimizing property of the ground state, i.e., that it satisfies a 
variational equation. This is done in Section 7 to estimate the L 1 -norm of the 
density in a region of the form {x : \x\ > R}. 

In Section 8 we establish the semiclassical estimates that allow one to 
compare the HF model with the TF model. To be more precise, there is no 
semiclassical parameter in our setup, but we derive bounds that in a semiclas- 
sical limit would be asymptotically exact. 

It turns out to be useful to use the electrostatic energy (or rather its square 
root) as a norm in which to control the difference between the densities in TF 
and HF theory. The properties of this norm, which we call the Coulomb norm, 
are discussed in Section 9. Sections 4-9 can be read almost independently. 

In Section 10 we state and prove the main technical tool in the work. It is 
a comparison of the screened nuclear potentials in HF and TF theory. Using a 
comparison between the screened nuclear potentials at radius R one may use 
the result of the separation of the outside from the inside given in Section 6 to 
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get good control on the outside region {x : \x\ > R}. Using an iterative scheme 
one establishes the main estimate for all R. The two main technical lemmas 
are proved in Section 11 and Section 12 respectively. 

Finally the main theorems are proved in Section 13. 

The main results of this paper were announced in [26] and a sketch of the 
proof was given there. The reader may find it useful to read this sketch as a 
summary of the proof. 

2. Notational conventions and basic prerequisites 

We shall throughout the paper use the definitions 

(11) B(r) := {y G M 3 : \y\ < r}, 

(12) B(x,r) := {y G R 3 : \y - x\ < r}, 

(13) A(n,r 2 ) := {x G R 3 : r x < \x\ < r 2 }. 

For any r > we shall denote by X r the characteristic function of the ball 
B(r) and by Xr = 1 — Xr- We shall as in the introduction use the notation 
[t]± = (t)± := max{±t,0}. 

Our convention for the Fourier transform is 

(14) f(p):=(2Try 3 / 2 je^f(x)dx. 
Then 

(15) fTg = {2vf 2 fg, ||/|| 2 = ||/|| 2 , \f{p)\<(2^/ 2 \\f\\ l 
and 

(16) J f f(x)\x-y\- 1 gjy)dxdy = 2(27i) j f(p)W)W 2 dp. 

Definition 2.1. (Density matrix). Here we shall use the definition that a 
density matrix, on a Hilbert space TC, is a positive trace class operator satisfying 
the operator inequality < 7 < /. When H is either L 2 (M 3 ) or L 2 (R 3 ; C 2 ) we 
write 7(2, y) for the integral kernel for 7. It is 2 x 2 matrix valued in the case 
L 2 (M 3 ; C 2 ). We define the density < p 7 G L 1 (M 3 ) corresponding to 7 by 

(17) /> 7 :=]>»,(x)| 2 , 

j 

where Vj and Uj are the eigenvalues and corresponding eigenfunctions of 7. 
Then j* p 7 = Tr[ 7 ]. 

Remark 2.2. Whenever 7 is a density matrix with eigenfunctions Uj and 
corresponding eigenvalues vj on either L 2 (M 3 ) or L 2 (M 3 ;C 2 ) we shall write 

(18) Tr[-A 7 ]:=]>> / \V Uj (x)\ 2 dx. 
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If we allow the value +00 then the right side is defined for all density matrices. 
The expression — A7 may of course define a trace class operator for some 7, 
i.e., if the eigenfunctions Uj are in the Sobolev space H 2 and the right side 
above is finite. In this case the left side is well defined and is equal to the right 
side. On the other hand, the right side may be finite even though — A7 does 
not even define a bounded operator, i.e., if an eigenfunction is in H 1 , but not 
in H 2 . Then the sum on the right is really 



Tr 



(_A) 1 /2 7 (_ A )i/2l =T r[V- 7 V]. 



It is therefore easy to see that (18) holds not only for the spectral decompo- 
sition, but more generally, whenever 7 can be written as 7/ = J2j u j( u ji f) u j: 
with < Uj (the Uj need not be orthonormal) . The same is also true for the 
expression (17) for the density. 

Proposition 2.3 (The radius of an infinite neutral HF atom). The map 
7 1 ^ Tr[— A7] as defined above on all density matrices is affine and weakly 
lower semicontinuous. 

Proof. Choose a basis /i,/2, ■ ■ ■ for L 2 consisting of functions from H 1 . 
Then 

Tr[-A 7 ] = ^(V/ m , 7 V/ m ). 

m 

The affinity is trivial and the lower semicontinuity follows from Fatou's lemma. 

□ 

We are of course abusing notation when we define Tr[— A7] for all density 
matrices. This is, however, very convenient and should hopefully not cause 
any confusion. 

If V is a positive measurable function, we always identify V with a mul- 
tiplication operator on L 2 . If Vp-y G L 1 (M 3 ) we abuse notation and write 

V 



7 

Tr[y 7 ] := J Vf> 



As before if happens to be trace class then the left side is well defined 
and finite and is equal to the right side. Otherwise, we really have / Vp^ = 



Tr 



lV2 r T/1 l/2 



Lemma 2.4 (The IMS formulas). If u is in the Sobolev space H l {R 3 ; C 2 ) 
or H 1 (M. 3 ) and if S G C 1 (1R 3 ) is real, bounded, and has bounded derivative 
then 1 



(19) Re J V (~V) • Vu = J \V(Eu)\ 2 - J |VS| 



2 1 |2 
It . 



1- We denote by u* the complex conjugate of u. In the case when u takes values in C 2 this refers 
to the complex conjugate matrix. 
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If 7 is a density matrix on L 2 (R 3 ; C 2 ) or L 2 (R 3 ) and if Si, . . . , E m G C 1 ( 
are reaZ, bounded, have bounded derivatives, and satisfy S 2 + . . . + z? m = 1 then 



(20) 



Tr [-A 7 ] 



Tr[-A(S l7 S 1 )]-Tr (V~!) 2 7 



+ ... 



+ Tr [-A(S m7 S m )] - Tr [(VS m ) 2 7 



Note that Sj 7 Sj again defines a density matrix (where we identified Sj with a 
multiplication operator) . 



Proof. The identity (19) follows from a simple computation. If we sum 
this identity and use S 2 + . . . + S 2 ^ = 1 we obtain 

J |Vu| 2 = J |V(Sm)| 2 -| |VS 1 | 2 M 2 + ... + J |V(S m n)| 2 - J \VE m \ 2 \u\ 2 . 

If we allow the value +oo this identity holds for all functions u in I? . Thus 
(20) is a simple consequence of the definition (18). □ 



Theorem 2.5 (Lieb-Thirring inequality). We have the Lieb- Thirring 
inequality 



(21) 



Tr 



±A 7 1 > K x \ p. 



9 5/3 



where K\ := 20.49. Equivalently, If V £ L 5 / 2 (R 3 ) and if 7 is any density 
matrix such that Tr[— A 7 ] < oo we have 



(22) 



where L\ := \ 



Tr 

\ 2/3 



-iA 7 ]-Tr[y 7 ]>-L!|[y]f , 



0.038. 



The original proofs of these inequalities can be found in [18]. The con- 
stants here are taken from [7]. From the min-max principle it is clear that the 
right side of (22) is in fact a lower bound on the sum of the negative eigenvalues 
of the operator — |A — V. 



Theorem 2.6 (Cwikel-Lieb-Rozenblum inequality). If V e L 3 / 2 (R 3 ) 
then the number of nonpositive eigenvalues of — ^A — V, i.e., 



Tr 



X(-oo,0] 



where X(-oo,o] * s the characteristic function of the interval (— oo,0], satisfies 
the bound 



(23) 



Tr 



X(-oo,o] (-h A ~ V ) 



[V] 



3/2 
+ ' 



where L := 2 3 / 2 0.1156 = 0.3270. 
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The original (independent) proofs can be found in Cwikel [4], Rozen- 
blum [19], and Lieb [9]. The constant is from Lieb [9]. 



3. Hartree-Fock theory 



In Hartree-Fock theory, as opposed to Schrodinger theory, one does not 
consider the full TV-body Hilbert space /\ N L 2 (R 3 ;C 2 ). One rather restricts 
attention to the pure wedge products (Slater determinants) 

(24) ^ = (iV!)- 1 / 2 ^ A ... Au N , 

where ui,...,un G L 2 (M. 3 ; C 2 ). Then one minimizes the energy expectation 

(*,tfjv,z*) 



of the Hamiltonian 

N 

(25) H NtZ := £ 

i=i 



oh E 



l<i<j<iV |Xi X ^ 



over wave functions * of the form (24) only. 

If 7 is the projection onto the iV-dimensional space spanned by the func- 
tions u±, . . . , ujy, the energy depends only on 7. In fact, 



= ^ HF (7). 



Here we have defined the Hartree-Fock energy functional 



(26) £ HF ( 7 ) : 



Tr 
Tr 



7 



+ P( 7 ) - £*( 7 ) 



-§A 7 ] ~y Z|x| _1 p 7 (s)<ic + I>(7)-£*(7), 

where we have introduced the direct Coulomb energy, defined in terms of the 
Coulomb inner product D (see also (79) below), by 

(27) D( 7 ) :=D(p 7 ,p 7 ) = 5// P 7 (x)|:r-y|- 1 p 7 (y)cZxd 2 / 

and the exchange Coulomb energy 



(28) 



£*( 7 ) : =\J J Tr c 2 l7(^,y)| 2 ^xdy. 



Definition 3.1. (The Hartree-Fock ground state). Let Z > be a real 
number and iV > be an integer. The Hartree-Fock ground state energy is 

E™(N,Z) :=inf{£: HP (7):7*=7, 7 = 7 2 , Tr[ 7 ] = n}. 

If a minimizer 7 HF exists we say that the atom has an HF ground state described 
by 7 HF . In particular, its density is P HF (x) = Pjhf(x). 
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Theorem 3.2 (Bound on the Hartree-Fock energy). For Z > and any 

integer N > we have 

E HF {N,Z) > -3(4 7 rL 1 ) 2 / 3 Z 2 iV 1 / 3 , 
where L\ is the constant in the Lieb- Thirring inequality (22). 

Proof. Let 7 be an N dimensional projection. Since the last term in H^^z 
is positive we see that £ HF (7) > Tr ^— — 7 . It the follows from 

the Lieb-Thirring inequality (22) that for all R > we have 

£ HF (7) > -Li [ Z 5 / 2 |x|" 5 / 2 dx - ZNR- 1 . 

J\x\<R 

The estimate in the theorem follows by evaluating the integral and choosing 
the optimal value for R. □ 

Remark 3.3. The function N h-> E nF {N,Z) is nonincreasing. This can 
be seen fairly easily by constructing a trial iV + 1-dimensional projection from 
any iV-dimensional projection by adding an extra dimension corresponding to 
a function u concentrated far from the origin and with very small kinetic energy 
/ |Vu| 2 . This trial projection can be constructed such that it has an energy 
arbitrarily close to the original iV-dimensional projection. Therefore we also 
have that 

(29) E HF (N, Z) = inf {s HF (~f) : 7* = 7, 7 2 = 7, Tr 7 < TV"}. 

This Hartree-Fock minimization problem was studied by Lieb and Simon 
in [16]. They proved the following about the existence of minimizers. 

Theorem 3.4 (Existence of HF minimizers). If N is a positive integer 
such that N < Z + 1 then there exists an N- dimensional projection 7 HF mini- 
mizing the functional £ HF in (26), i.e., E HF (N,Z) = £ HF (7 HF ) is a minimum. 



In the opposite direction the following result was proved by Lieb [13]. 

Theorem 3.5 (Lieb's bound on the maximal ionization). If N is a 

positive integer such that N > 2Z + 1 there are no minimizers for the Hartree- 
Fock functional among N-dimensional projections, i.e., there does not exist an 
N- dimensional projection 7 such that <? HF (7) = E HF (N, Z). 

This theorem will, in fact, follow from the proof of Lemma 7.1 below (see 
page 503). Although this result is very good for Z = 1 it is far from optimal 
for large Z. In particular the factor 2 should rather be 1. This fact known as 
the ionization conjecture is one of the of the main results of the present work. 
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Theorem 3.6 (Universal bound on the maximal ionization charge). There 
exists a universal constant Q > such that for all positive integers satisfying 
N > Z + Q there are no minimizers for the Hartree-Fock functional among 
N- dimensional projections. 

Remark 3.7. Although, it is possible to calculate an exact value for the 
constant Q above it is quite tedious to do so. Moreover, the present work 
does not attempt to optimize this constant. The result of this work is mainly 
to establish that such a finite constant exists. This of course raises the very 
interesting question of finding a good estimate on the constant, but we shall 
not address this here. 

The proof of Theorem 3.6 is given in Section 13 on page 534. 

Theorem 3.8 (Bound on the ionization energy). The ionization energy 
of a neutral atom E HF (Z—1, Z) — E HF (Z, Z) is bounded by a universal constant 
(in particular, independent of Z). 

This theorem is proved in Section 13 on page 573. 

The variational equations (Euler-Lagrange equations) for the minimizer 
was also given in [16]. Since the Hartree-Fock variational equations shall be 
used later in this work, we shall derive them in Theorem 3.11 below. 

We first note that the Hartree-Fock functional £ HF may be extended from 
projections (i.e., density matrices with 7 2 = 7) to all density matrices. If 
Tr [— A7] < 00 all the terms of £ HF are finite. In fact, Tr [Z|x| _1 7] is finite 
by the Lieb-Thirring inequality (21) since Z\x\~ x G L°°(IR 3 ) + L 5 / 2 (M 3 ). The 
term T>("f) is finite by the Hardy-Littlewood-Sobolev inequality since p 7 £ 
L X (IR 3 ) n L 5 / 3 (M 3 ) C L 6 / 5 (M 3 ). Finally, £x(j) < D(j) since 

Din) - &(7) = 3 E / / 1 _ 1 W " C ® C dy > 0, 

J J \x y\ 

1,3 

when Ui are the eigenvalues of 7 with U{ being the corresponding eigenfunctions. 
If Tr [— A7] = 00 we set £ HF (7) := 00. It is clear that lim n £ HF (7 n ) = 00 if 
lim n Tr [-A7J -> 00. 

Remark 3.9. It is important to realize that although .0(7) — £x(~j) is 
positive it is not a convex functional on the set of density matrices. In partic- 
ular, the Hartree-Fock minimizer need not be unique. (A simple example of 
nonuniqueness occurs for the case N = 1. For a one-dimensional projection 7, 
it is clear that D(^) — £x(pj) = 0, hence the minimizer in this case is simply 
the projection onto a ground state of the operator — |A — Z\x\^ on the space 
L 2 (M 3 ;C 2 ). There are many ground states since the spin can point in any 
direction.) 
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Another fact related to the nonconvexity of the Hartree-Fock functional 
is the important observation first made by Lieb in [11] that the infimum of 
the Hartree-Fock functional is not lowered by extending the functional to all 
density matrices. For a simple proof of this see [1]. 

Theorem 3.10 (Lieb's variational principle). For all nonnegative inte- 
gers N we have 

inf{£ HP ( 7 ): 7 *= 7 , 7 = 7 2 , Tr[ 7 ] = n} 

= inf {£ HF ( 7 ) : < 7 < /, Tr[ 7 ] = N} 

and if the infimum over all density matrices (the inf on the right) is attained 
then so is the infimum over projections (the inf on the left). 

We now come to the properties of the Hartree-Fock minimizers, especially 
that they satisfy the Hartree-Fock equations. These equations state that a min- 
imizing ^-dimensional projection 7 HF is the projection onto the iV-dimensional 
space spanned by eigenfunctions with lowest possible eigenvalues for the HF 
mean field operator 

(30) F 7 hf := -±A - Z\x\~ l + p UF * \x\- 1 - /C 7 hf. 

Here /C 7 hf is the exchange operator defined by having the 2 x 2-matrix valued 
integral kernel 

/C 7 H F (x,y) := \x - y\~ 1 7 UF (x,y). 

Thus 7 HP (x,y) = Y,iLi u i( x ) u i(y)*, where H^fUi = em, and E\, e 2 , . . . , £jv 
< are the N lowest eigenvalues of H^rf counted with multiplicities. 

This self-consistent property of a minimizer 7 HF may equivalently be stated 
as in the theorem below. 

Theorem 3.11 (Properties of HF minimizers). 7/ 7 HF with density p HF is a 
projection minimizing the HF functional £ HF under the constraint Tr [ 7 HF ] = N 
then p HF G L 5 / 3 (1R 3 ) n L 1 (IR 3 ) and H^kf defines a semibounded self -adjoint 
operator with form domain // 1 (IR 3 ;C 2 ) having at least N nonpositive eigen- 
values. Moreover, 7 HF is the N- dimensional projection minimizing the map 
7 ^ Tr H^hf'j . 

Remark 3.12. The reader may worry that, because of degenerate eigen- 
values of H^hf, the iV-dimensional projection 7 minimizing Tr H^hf^ may 
not be unique. That it is, indeed, unique was proved in [2]. 

Proof of Theorem 3.11. We note that Tr [7 HF ] = N, Tr [-A7 HF ] < 00, 
and the Lieb-Thirring inequality (21) implies that p HF G L 5/3 (M 3 ) n L 1 (R 3 ). 
From this it is easy to see that p HF * is a bounded function (in fact, it 
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is continuous and tends to as \x\ — > oo). Moreover, in the operator sense 
/C 7 hf < p HF * |x| _1 . This follows, since for / G L 2 (M 3 ;C 2 ) we have 

/ P m *\x\- l \f{x)\ 2 dx- J J f(x)*K^(x,y)f(y)dxdy 



N 



\ui(x) <8> f(y) - f(x)0Uj(y)\\ 
\x-y\ 



-dx dy, 



where u\,...,un is a complete set of eigenfunctions of 7 HF . It is therefore 
clear that H^hf defines a semibounded operator with form domain if 1 (IR 3 ; C 2 ). 



Thus it makes sense to compute Tr 



H, 



7 



hf7 



if and only if Tr [— A7] < 00. 



Let 7' be an iV-dimensional projection with Tr [— A7 ; ] < 00. We shall 
prove that 



Tr 



> Tr 



# 7 hf7 hf 



For < t < 1, consider the density matrix 7$ = (1 — i)7 HF + £7'. It 
satisfies Trfyt] = N. By the Lieb variational principle, Theorem 3.10, we have 
that £ HF (7 HF ) = £ HF (7o) < £ UF (jt), for all < t < 1. Hence 

^ - — — = Tr //- 1:1 - — Tr 



dt 



t=o 



i? 7 HF7 HF 



The fact that Tr 



# 7 hf7 



is minimized among TV-dimensional projections 



implies in particular that H^hf has at least N nonpositive eigenvalues. □ 



4. Thomas-Fermi theory 

In this section we discuss the facts needed from Thomas- Fermi theory. We 
focus only on the results that we shall use in our study of Hartree-Fock theory. 

Definition 4.1. (Thomas-Fermi functional). Let V G L 5 / 2 (M 3 ) + L°°(1R 3 ) 

with 

inf {llWHioofRS) : V — W £ L 5/2 (M 3 )} = 0. 
Corresponding to V we define the Thomas-Fermi (TF) energy functional 

4 f (p) = U^ 2 f" j p 5/3 -f v P + ^fl p^\ x ~ y\~ 1 p(y)dxd y , 

on functions p with 0<pG L 5 / 3 (M 3 ) n L 1 (R 3 ). 

Note that the Hardy-Littlewood-Sobolev inequality implies that D(p, p) = 
\ I I p(x)\x — y\~ l p(y)dx dy is finite for functions p G L 5 ^ 3 DL 1 C L 6 / 5 . Hence 

TF 

£,/ is finite on these functions. 
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multiplier) p^(N'), with < p^(N') < supV, such that p^J is uniquely 



9 tf 
v 



The proof of existence and uniqueness of minimizers to the TF functional 
and the characterization of their properties can be found in the work of Lieb 
and Simon [17] (see also [10]). We state the properties that we need in the 
following theorem. 

Theorem 4.2 (The TF minimizer). Let V be as in Definition 4.1. For all 
N' > there exists a unique nonnegative G L 5 / 3 (M 3 ) such that J p^J < N' 
and 

(31) £y{p™) = inf |4 F (p) : p G L 5 / 3 (M 3 ), j p < AT' J. 

On the other hand there exists a (unique) chemical potential (Lagrange 
tiplier) A t y F ( J 
characterized by 

(32) £ T v Y (p^) + ^(N')f p\ 

= inf \£ T y(p) + p T v F (N') J p : < p G L 5 / 3 (M 3 ) D ^(M 3 ) 

Moreover, p^J is the unique solution in L 5 / 3 n L 1 to the Thomas-Fermi 
equation (the Euler- Lagrange equation for the variational problem (32)) 

(33) ±(3vr 2 ) 2 / 3 (p-(x)) 2/3 = [V(x) - p™ * M" 1 - p^(N')] + . 

Ifp TF (N') > thenjp T v F = N f . Therefore p TF (N')fp TF = p TF (N')N'. 

For all < p there is a unique minimizer < p G L 5 / 3 n L 1 to £^ F (/9) + 
jtx/p. (If p > supV then p is simply zero.) 

We shall be interested in properties of the Thomas-Fermi potential 

(34) ip™ := V(x) - p™ * \x\-\ 

The Thomas-Fermi equation (33) can be turned into the Thomas-Fermi dif- 
ferential equation 

(35) = 2^(3^1 \<p™ - ^(N')} 3/2 + AV, 



v 

which holds in distribution sense. 



Theorem 4.3 (Maximal ionization). There exists a nonnegative real 
number N c , possibly equal to +oo, such that p^(N') > if and only if N' < 
N c . Moreover, 

(36) N c > liminf(4vr)" 1 / rV(ru)duj, 

r^oo J §2 

where duo is the surface measure on the unit 2- sphere S 2 . 
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Proof. Since £ v is a convex functional of p it is clear that £\, 



P 



is 



a convex and decreasing function of N'. Hence there is a value N c such that 
£y F (Py F ) is strictly decreasing for N' < N c and constant for N' > N c . Thus 
if N' > N c then / p™ = N c . Since / p™ = N' if p^J > we must have 
H™(N') = p^(N c ) = for N' > N c . On the other hand since f = N' 
if N' < N c we cannot have p^(N') = in this case. This proves the first 
assertion. 

In order to prove the second assertion we may of course assume that 
N c < oo. Since p^(N c ) = we have for the corresponding Thomas-Fermi 
minimizer that 



/ p T J{ru)du = (Const.) / \v(ru) - p 
Js 2 v Js 2 L 



5™ * |roj| _1 



3/2 
+ 



duj 



> (Const.) 



Utt)- 1 [ V(ru)duj -r' 1 ( p 
Js 2 Jr 3 



3/2 
+ 



where the last estimate follows from Jensen's inequality and Newton's theorem. 
Since we are considering a TF minimizer p^J such that / p^J = N c it is clear 

that if (36) is violated then J s2 p^J (ru)duj > cr~ 3 / 2 for some positive constant c 
and all large enough r. Hence N c = J p^J = oo in contradiction with our 
assumption. □ 

Proving a bound on N c in the opposite direction is in general more difficult. 
We shall return to a partial converse to (36) in Corollary 4.8 below. 

Usually the Thomas-Fermi model is studied for the potential V being 
the Coulomb potential, i.e., In this case we denote p^J , <p^ F , and 

p^J simply by p TF , (p TF , and p TF . These are the functions discussed in the 
introduction. In fact, the equations (5) and (6) correspond to (33) and (35). 

From Theorem 4.3 we see that in this case N c > Z. We shall see below 
after Corollary 4.8 that indeed N c = Z. 

The first mathematical study of the atomic TF equation was done by 
Hille [6]; a much more complete analysis can be found in [17]. 

The function ip TF satisfies the asymptotics ip TF {x) ~ 3 4 2 -3 7r 2 |x| for 
large x. The important thing to note about this asymptotics, first discovered by 
Sommerfeld [27], is that it is independent of Z. The Sommerfeld asymptotics 
is central to the present work and we shall prove a strong version of it in 
Theorems 5.2 and 5.4 below. Similar asymptotic estimates may be derived 
for the density using the TF equation (33). We shall more generally prove 
asymptotic bounds for ip^ , in the case when the potential V is harmonic in 
certain regions of space. 
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We now come to the main technical lemma in this section, which is a 
version of the Sommerfeld estimate. 2 

Lemma 4.4 (Sommerfeld estimate). Assume that ip > is a smooth 
function on \x\ > R and satisfies the differential equation 

Aip(x) = 2 7/2 (3Ti)- 1 (p(x) 3/2 , for \x\ > R, 

for some R>0. Let ( := (-7 + y/73)/2 « 0.77. Define 

(p(x) > 1/2 



and 



a(R) := liminf sup 

r\R \ x \ =r 



A(R) := liminf sup 

r\R | x | =r 



3 4 2 -3 7r 2 r -4 

3 4 2- 3 vr 2 r- 4 " 



- 1 



Then for \x\ > R we have 

(37) (l + a(R)\x\-<y 2 < 



<p(x) 



3 4 2- 3 7T 2 |x|- 4 



< 



(l + A(R)\x\-<) . 



Remark 4.5. It is important to realize that we are not assuming that ip 
is spherically symmetric. The lemma above can therefore not be proved by 
ODE techniques. By elliptic regularity the smoothness of ip would of course 
be a consequence of a much weaker assumption. 



Proof of Lemma 4.4. We first prove that <p(x) — > as \x\ — > oo. For 
this purpose consider L > 4R and for L/4 < \x\ < L the function f(x) = 
C[(\x\ - L/4)- 4 + (L - |x|)" 4 ]. We compute 



A/ 



C 



20 



x 



L/4)- e + (L-\x\)- e ) 



+ 8\x\~\L - \x\y 5 - 8\x\-\\x\ - L/4)- 5 
< 44C(L-\x\)- 6 + 20C(\x\- L/4)' 6 . 



On the other hand, f{xf/ 2 > C 3 / 2 ((\x\ - L/4)- 6 + (L - |x|)- 6 ). It is 
therefore clear that we can choose C independently of L such that A/ < 
2 7 / 2 (3vr)" 1 / 3 / 2 . We claim that tp(x) < f(x) for L/4 < \x\ < L. This is trivial 
for | a; | close to L/4 or close to L since here f(x) diverges whereas ip(x) remains 



2 A version of this Sommerfeld estimate was stated in the announcement [26]. The result stated 
was weaker than here in the sense that the exponents in the error terms were different for the upper 
and lower bounds. The result in the announcement also contained a minor error because the lower 
bound had been stated incorrectly. The better and correct version is the one stated and proved here. 
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bounded. Consider the set {L/4 < \x\ < L : ip(x) > f(x)}. This is an open set 
on which A((p - /) > 2 7 / 2 (37r)~ 1 (^ 3 / 2 - f 3 / 2 ) > 0; i.e., <p - f is sub-harmonic 
on the set and is zero on its boundary. Hence <f(x) < f(x) on the set which 
is a contradiction unless the set is empty. Thus for all L > 4R we have 
sup N=L/2 <p(x) < C ((1/4)" 4 + (1/2)- 4 ) L -4 . Hence, ^(x^x) 4 is bounded. 

Next we turn to the proof of the main estimate. Let R' > R and set 
A' = A(R') and a' = a(R'). Then a' and A' are finite. We consider the two 
functions 

U3%{x) := 3 4 2~V|x|- 4 (l + A'\x\- ( ) 

and 

u-,{x) := 3 4 2-V|x|- 4 (l + a'|x|-^)- 2 . 

Note that by the definition of a' and A' both functions are well-defined and 
positive for \x\ > R' . We claim that 

(38) Au+, (x) < 2 7 / 2 {3tt)- 1 uj+ {xf 2 and Au~,{x) > 2 7 / 2 {3tt)- 1 lj^{x) 3 / 2 . 

As we shall first show the lemma is a simple consequence of the estimates 
in (38). We give the proof for the upper bound. The lower bound is similar. 
Let 

Q + := > R' : ip(x) > u/£,(x)j. 

On Q + , if — uj\, is subharmonic. On the boundary of f2+, if — u>\, vanishes. 
For the subset n {x : |x| = R'} this follows from the choice of A'. Since 
(p(x) and uj\, (x) both tend to zero as \x\ tends to infinity we conclude that 

n + = $. 

Therefore ip(x) < oj~^, r ,Jx) for |x| > R' . For \x\ > R we get ip(x) < 
liminf R / Xi? w+ (R/) (x) = u\ {R) {x). 

It remains to check (38). For uj~, we get 

A^,(x) = 2 7 / 2 (3^)- 1 cu a ;(x) 3 / 2 ( l + (l- I C (C + 7))a'|x|^ 

+i(l + a'|x|-0- 1 (C«>l" C ) 2 )- 

Since C(C+ 7 ) = 6 and l+a'\x\-^ > we see that Aoj~, (x) > 2 7 / 2 (3vr)~ 1 w a ;(x) 3 / 2 . 
For we have 

Au+{x) 

= 2 7 / 2 (3vr)-^+ (x) 3 / 2 (l + A'\x\-t)- 3 / 2 (l + (l + ^J^) A'\x\~^ 

<2 7 / 2 (3tt)- 1 o;+(x) 3 / 2 , 
where we have used that 

(1 + A'\x\-^ 2 > 1 + lA!\x\~t = ! + (! + ^C(C + 7))A'\x\-<. □ 
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We can immediately use this lemma to get estimates on ip^ when //^ F = 0. 
For general jjl^ the result can be generalized as follows. 

Theorem 4.6 (Sommerfeld estimate for general l^)- Assume that V 
is continuous and harmonic for \x\ > R and satisfies lim| x |_ >00 V(x) = 0. 
Consider the corresponding Thomas- Fermi potential ip^ , which satisfies the 

TF differential equation (35). Assume that ^ < liminf inf ipJJ{x). Define 



r\R \x\=r V 



(39) 

and 
(40) 



a(R) 



liminf sup 

r\R | x | =r 



TF/ \ \ -1/2 
<Py (X) \ 



A(R, n 



TF^i 

V > 



liminf sup 

r\R | x | =r 



3 4 2 -3 7r 2 r -4 



3 4 2- 3 vr 2 r" 4 

T/ien acam, urctfi C = (-7 + \/73)/2 « 0.77, we /md /or a// \x\ > R 

(41) 
and 



■ iF >r c )+/4 F 



(42) 

where 

(43) 



<^( x ) < 3 4 2- 3 vr 2 |x|- 4 (l + A(i?, /i 
^ F (x) > max|3 4 2- 3 7r 2 |xr 4 (l + a(i?)|xr c )^ 2 , 
i/(/^ F ) := inf fl max|3 4 2- 3 7r 2 |xr 3 (l + a(R)\x\^y 2 ,^ F |x|| 



Proof. Since G L 5 / 3 PlL 1 it is easy to see that p^ F * |x| _1 is continuous 
and tends to zero as x tends to infinity. Thus from the assumption on V it 
follows that ip^J is continuous on \x\ > R and satisfies <p^{x) — ► as \x\ — > oo. 

Let i?' > i? and set A' = A(i?', ^ F ) and a' = a(R'). Then a' is well-defined 
if R' is close enough to R since then we may assume that tp™{x) > [i™ > 
for all |x| = R' . Both a' and A' are finite. Using the notation from the proof 
of Lemma 4.4 we define 

t (x) := uj\, (x) + and w^ TF ; (x) := max joT",(x), ^'|x| _1 | , 



Ct> TF 

where 

.' 

\x\>R' 



v' :— >i:iii max j|x|o; a ,(x), |x|/^ F j . 



Note that, since we assume that <p^(x) > ji^J for |x| = i?', we have that both 



y \ I t~y 

(x) and co~,(x) are positive for all \x\ > R'. We also have that w~i(x) > fiy F 



for \x\ = R' and hence that oj a ,(xo) = n^J at points xo where the minimum, 
defining u', is attained. (Note that |x|u/^(x) is a radially decreasing function 
for \x\ > R'.) 
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The proof of the present lemma is now similar to that of Lemma 4.4 if we 
can show that for Ixl > R' 



(44) 
and 
(45) 



Aw+ 4/ (x)<2 7 /2(3vr)- 1 

V 



Auj TF 

My 



"„TF (X) 



/i 



V 



> 2 7 / 2 (3vr)- 1 



1 3/2 



J + 



3/2 



Mt/ > 



M, 



(in distribution sense). The inequality (44) follows immediately from the first 
inequality in (38). The inequality (45) is slightly more complicated. Note 

v'\x 



that the definitions of lo TF and of u ! imply that to T f (x 

My My 



if 



U TF < V 

My > a ' 



and uj~ 



0J a ,{x) if w tf ,( x ) > V^v ■ Thus if 
u;~ T f ,(#) < My F we have that w _ TP f is harmonic. Hence (45) holds in 



My - ' ■ ^ 

this region. If go 



My 



(:r) > /i^ F then w TP (x) = oj n ,(x) and (45) follows 

My > a ' ^ My : a 

in this region from the second inequality in (38). Finally, since the maxi- 
mum of two subharmonic functions is also subharmonic, it is clear that the 
distribution Aw" TP is a positive measure and in particular positive on the 

My .a' 

set (of Lebesgue measure) zero where oj tf ,(^) — ~ Hence, (45) holds in 

My - a ' 



distribution sense for all \x\ > R'. 



□ 



As an application of the lower bound on ip^J in (42) we can get an estimate 



v 

on the chemical potential ^ . 

Corollary 4.7 (Chemical potential estimate) 



With the assumptions 

TF 

V 



and definitions in Theorem 4.6, in particular, if < liminf r ^flinf| a; | =J , <p^ F (x) 
we have 
(46) 



o3/4 / , 

(My F ) 3/4 < —^{l+\a{R)\R- C f /2 (j^r 1 j s2 rV{ru)du 



Pl F (y)dy 



Proof. According to (42) we have v{[iy F ) < liminfj, r |^ 00 |x|(/3^. F (x). Using 
that V is harmonic and tends to zero at infinity we have that for all r > R 

liminf \x\^(x) < (47T)- 1 / rV{ruj)duj - ( p T J{y)dy. 
\x\->oo v Js 2 JR 3 V 

Moreover since, ^ > the assumption /x^. F < lim inf j-x^ inf | x | =r (f v F (x) im- 
plies that the spherical average of V is nonnegative. 
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On the other hand, since (l + \a(R)\R~^ 2 < (l + a(R)\x\~^ 2 for 
\x\ > R, we have from (43) that K/^y F ) > v', where 

v 1 = minmax|3 4 2- 3 ^ 2 |xr 3 (l + |a(-R)l^) ^/^m} 

= 3-2- 3 /V/ 2 (l + |a(i?)|i?^)" 1/2 (4 F ) 3 / 4 . □ 

This corollary immediately gives a partial converse to Theorem 4.3. 

Corollary 4.8 (Upper bound on maximal ionization). IfVis harmonic 
and continuous for \x\ > R and satisfies V(x) — > as \x\ — > oo and if moreover 
< limmi r \ R mi\ x \ =r ip™ (x) then 

I P T J < lim Utt)- 1 f rV(ru:)duj. 

J V r->oo J S 2 

In particular, if liminf r \ N /jinf| :r | =r (p^J(x) > (which may not necessarily be 
true) we have 

N c < Inn (47T)" 1 f rV(ru)dw. 



Remark 4.9. The limit above of course exists since by the harmonic- 
ity of V and since V tends to zero at infinity we have that f S 2 rV(ruj)duj is 
independent of r. 

The difficulty in using Corollaries 4.7 and 4.8 in concrete examples lies in 
establishing the condition 




5. Estimates on the standard atomic TF theory 

In the usual atomic case the Coulomb potential V(x) = Z\x\~ x is harmonic 
away from x = and we can use Corollary 4.8 for all R > 0. Since p TF * 
is a bounded function it follows that (p TF (x) — > oo as x — > 0. The condition 
(47) is therefore satisfied if R is chosen small enough. It therefore follows 
from Theorem 4.3 and Corollary 4.8 that N c = Z. Thus the neutral atom 
corresponds to p7 F = 0. 

LEMMA 5.1. Let (f^ F be the TF potential for the neutral atom then if <fi TF 
is the potential for a general /x TF > we have 

fo (x) <<p {x)< <p {x) + n . 
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Proof. See Corollary 3.8 (i) and (iii) in [10]. □ 

We now easily get an upper bound agreeing with the atomic Sommerfeld 
asymptotics. 

Theorem 5.2 (Atomic Sommerfeld upper bound). The atomic TF po- 
tential satisfies the bound 

<p TF (x) < min{3 4 2- 3 7r 2 |xr 4 + p TF , Z^ 1 }. 



Proof. This follows immediately from and (34) and (41) together with the 
fact that p TF is nonnegative. Simply note that since v? TF (x)|x| — > Z as x — > 
we have that A(0, p TF ) = in (41). □ 

Lemma 5.3 (Lower bound on the TF potential). In the atomic case we 
have for all N > and Z > 

^ F (x) > Zlxr 1 - min {N\x\-\ J^^ 1 *} ■ 



Proof. We have by Newton's theorem 

p™*^ 1 = \x\- 1 [ P TF (y)dy+( P^iyM^dy 

J\y\<\x\ J\y\>\x\ 
< mmj^lxl- 1 , J P^iyM^dy} . 

From the Sommerfeld upper bound Theorem 5.2 and the TF equation (33) we 
have 

p TF (x) 2 / 3 < 2(3vr 2 )- 2 / 3 min{3 4 2-V|x|- 4 , Z^ 1 }. 

Hence 

p TF (x) <mm{ Cl Z^ 2 \x\- 3 / 2 , c 2 \x\- 6 } , 

where a := 2 3 / 2 (3vr 2 )- 1 and c 2 := 3 5 2- 3 vr. Let r := (c 2 / 'ci) 2 / 9 Z' 1 / 3 . When 
\x\ = ro the two functions, in the minimum above, are equal. Thus 



/ P^ivM^dy < 4vr Cl Z 3 / 2 f ° t-^dt + 4vrc 2 /°° t- 5 dt = ±^cl /9 cl /9 Z 4 / 3 
J JO Jro 3 

Z 4 3 . 



(9vr) 2 / 3 

The lemma follows from the definition (34) of the TF potential. □ 
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Theorem 5.4 (Atomic Sommerfeld Lower bound). The TF potential 
satisfies 

( Z\x\~ x - 22(9tt)- 2 / 3 Z 4 / 3 , if \x\ < P0Z- 1 / 3 



ip TF (x) > < 



2 ■ 4 



max! 3 4 2- 3 ^ 2 fl + aZ^/ 3 |x|~^ |.r 

1 V 1 ' ' tf\x\>P Z-^ 

(Z-N^lxr 1 ^ 

where (5q = ^ 9 ^ 4 / and ( = (—7 + \/73)/2 as in Theorem 4.6 and a = 43.7. 

Proof. Let R = (9vr) 2 / 3 Z" 1 / 3 /44. Note that for \x\ < R the bound we 
want to prove is identical to the bound in Lemma. 5.3. 

If iV > Z, i.e., /U TF = the lower bound follows from Theorem 4.6 since 
a is chosen so as to make the lower bound continuous at \x\ = R and at these 
points we clearly have ip TF (x) > = /U TF . 

For general N the lower bound follows from the case N = Z because of 
Lemma 5.1 and Lemma 5.3. □ 

We end this section by giving a bound on the screened nuclear potential 
Q R at radius R in the atomic TF theory. 



Lemma 5.5 (Bound on $^ F ). We have 

■ TF 



<^>)< 3 4 2-vi*r 4 +^ 



Proof. We write ^JZ(x) = (p TF (x) + P TF (y)\x - y\ 1 dy. From 

Theorem 5.2 and the TF equation (33) we see that 

p^{y) = 2 3 / 2 (3^ 2 )- V F (y) - ^ TF ]f < 2- 3 3 5 7r|y|- 6 

and hence 

/ p™{y)\x-y\- l dy < f 2^^% - y^dy 

J\y\>\x\ J\y\>\x\ 

= [ 2- 3 3 5 7r|y|- 7 ^ = 2- 3 3V|x|- 4 . 

J\y\>\x\ 

The lemma follows from Theorem 5.2. □ 

6. Separating the outside from the inside 

We shall here control the energy coming from the regions far from the 
nucleus. Let 7 HF be an HF minimizer with Tr[7 HF ] = N. (We are thus assuming 
that iV is such that a minimizer exists.) 
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Definition 6.1. (The localization function). Fix < A < 1 and let 
G : R 3 -> R be given by 



G(x) 



if |x| < 1 

(ir/2)(\x\ - 1) [(1 - A)- 1 - I]" 1 if 1 < |x| < (1 - A)- 1 
k vr/2 if (1 - A)" 1 < \x\. 



We introduce the cut-off radius r > and define the outside localization 
function r {x) = sinG(|x|/r). Then 

= if \x\ < r 
0<e r (x)< <1 if r < \x\ < (1 - A) _1 r 
= 1 if (1 - A)-V < |x| 

and |W r (x)| 2 + |V(1 - 9 r (x) 2 y/ 2 \ 2 < (^/(2Ar)) 2 (since (1 - A)" 1 - 1 > A). 

We shall consider the HF minimizer restricted to the region {x : \x\ > r}. 
We therefore define the exterior part of the minimizer 



(48) 



e r 7 HF e r 



and its density p™(x) = 6 r (x) 2 p HF (x). In order to control 7™ we introduce 
an auxiliary functional defined on all density matrices with Tr [— A7] < 00 (see 
Remark 2.2) by 



(49) 



£: A ( 7 ) = Tr[(-iA-^ F ) 7 ]+i J J Pl (x)\x-y\- l Pl (y)dxdy, 



|a;|>T' 
\y\>r 



where the screened nuclear potential $ is defined in (3) in Definition 1.1. 
Note that the functional £ A , in contrast to the HF functional £ HF in (26), does 
not contain an exchange term. 

The main result in this section is that 7™ almost minimizes £ A . More 
precisely, we shall prove the following theorem. 



Theorem 6.2 (The outside energy). For all < A < 1 and all r > we 



have 



(50) £ A [7 r HF ] < inf jf A (7) : supp p 7 C {y : \y\ > r}, J Pj < J X^P HF } + U, 



where the error is 



(51) K = C x (r) 



p HF (x)dx + 2L 1 



$"l x)r (x)f 2 dx + £ x (X 



|x|>(l-A)r 



(l-A)r<|a;|<(l-A)- 1 r 
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with 



IT 2 7T 



Cx(r) = [ - , 2 + 



^(A(l-A)r) 2 r\J 
Here L\ is the constant in the Lieb- Thirring inequality (22). 

Proof. Besides 6 r we introduce two other localization functions 



(l-^(l-A)) 1/2 and 0(O) = (5 2 (1- 



A) ~ ®r 



1/2 



Note that 9 2 _ + 9 2 0) + 2 = 1 and that 

(W_) 2 + (V0 (o) ) 2 + (V9 r ) 2 < (vr/(2A(l - A)r)) 2 . 
We introduce the inside part of the HF minimizer 

(52) 7* F = 0-7 HF 0- 

We shall prove (50) by showing that for all density matrices 7 with supp p 1 C 
{x : \x\ > r} and / p 1 < J Xr P HF we have 

(53) £ A [7 r HF ] + £ HF [7 HF ] - n < £ HF [7 HF ] < £ A [l] + £ HF [7 ro ] , 
with 1Z given by (51). The estimate (50) follows immediately from (53). 

Proof of the upper bound in (53) . Since 7 HF is a minimizer for <? HF under 
the condition Tr[7 HF ] < N (see (29)) we have for any density matrix 7 with 
Tr [7] < N that £ HF (7 HF ) < £ HF (% We take 

Since the support of p 7 is disjoint from the support of 0- we see that 7^7 = 
and hence 7 is a density matrix. 
Note that 

Tr [7] = Tr [ 7 HF ] + j P 7 < Tr [ 7 ro ] + j xt P HF = J (e 2 _+xt)p HF < J P HF < N. 

We shall compute £ HF (7). The only terms in £ HF that are not linear in 
the density matrix (and thus do not simply split into a sum of terms for 7 HF 
and 7) are the exchange and direct Coulomb energies. Because of the support 

properties of 7 HF and 7 we have that 7 2 = (j HF ^j + 7 2 and therefore even the 
exchange term satisfies 

£ X (7) = £* (7^) + £x (7) > £x (7^ F ) • 
We are thus left with the direct Coulomb energy. For this we find that 

V (7) = V (7 HF ) + V (7) + j e 2 _{y) P « F (y)\x - y^P^dy. 
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By the choice of the support of we have that 

J 9 2 _{y) P ^(y)\x-y\- l P ^x)dxdy < J - <Z>™ (x)) P ^(x)d. 
We have thus proved the upper bound in (53). 



□ 



Proof of the lower bound in (53). Let again the inside part of the HF 
minimizer be 7^ F defined by (52) and introduce also the middle part 7 ( 
6 l (o)7 HF 6'(o)- Since 9 2 _ + 9 2 0) + 2 = 1 we have from the IMS formula (20) that 



(o) 



(54) 



Tr 



-±A 7 HF 



Tr 



i A ( 7 HF + 7nf + 7, HF 



(0) 



±Tr 

2 ±L 



r 



((w_) 2 + (w (0) ) 2 + (w r ) 2 ) 



> Tr 



-iA(7 HF + 7^ + 7, HF 



(tt 2 /8) (A(l - A)r)" 



p ilF (x)dx. 



(l-A)r<|x|<r(l-A)- 1 

We now come to the lower bounds on the Coulomb terms. Note that 



> 



(e 2 _{x) + e 2 0) (x) + e 2 (x)) (e 2 _{ y ) + e 2 0) ( y ) + e 2 ( y )] 

9 2 _(x)e 2 _(y) + 9 2 (x)9 2 (y) + 9 2 (x) (9 2 _{y) + # 2 } (y)) 
+ (9 2 _{x) + 9% } (x)) 9 2 (y) + 9 2 Q) {x) 2 9 2 _ (y) + 9 2 _(x) 2 9 2 a) {y) . 

Note that 9 2 _ (x) + 9 2 0) (x) > Xr(x) and 9 2 _ (x) > X(i-\) r We may therefore 
estimate the Coulomb kernel from below by 



x-y\ >V(x,y), 



where 



(55) V(x,y) := 9l (x)\x - y^ (y) + 9 2 (x)\x - y\-H 2 r (y) 
+9 2 (x)\x - y\- l Xr{y) + X r (x)\x - y\~ 1 9 2 r {y) 
+9%){x)\x - j/|" 1 X(i-A)r(j/) + X(i-x) r (x)\x - y\~ 1 9 2 0) (y). 

The function V is pointwise positive and symmetric in x and y. 

Recall that 7 HF is a projection onto the subspace spanned by the orthonor- 
mal vectors ui, U2, ■ ■ ■ , un and 
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(56) 



\ui(x) <g> Uj(y) - Uj(x) <g> Ui(y)\\ 
\x-y\ 



-dx dy 



P(7 HP )-^(7 HF ) = iE// 

= \ J J P HF (x)V(x,y) P 11F (y)dxd y 

Tr c2 \r F (x,y)\ 2 V(x,y)dxdy. 
We estimate these two terms independently. We obtain for the first term in (56) 
(57) \ J f p HF (x)V(x, y) P HF (y)dx dy = V ( 7 HF ) + V (7 r HF ) 

+/ (iff - ^ F(X) ) p * F{x)dx + Tt [(r - ^Um) ■ 

To estimate the last term in (56) we use that for all \x\ > r 

e r (xf < e r ( x ) < (7r/2)(|x|/r - 1) ((i - A)- 1 - . 

Thus if \y\ < r we have \x—y\~ 1 9 r (x) 2 < (7r/2)r _1 ((1 - A) -1 - l) -1 and hence 
J f Tr : , [| 7 HF (x,y)| 2 ] \x - y^ 1 (x r (y)0 2 r (x) + X (1 _ A)r (y)^ 0) (x)) dxdy 



< ^((l-A^-^^l + a-A)- 1 ) J J Tr c2 [| 7 HF (x,y)r 



dx dy 



\y\<r 
\x\>(l-X)r 



Tr c2 |7 HF (x,y)r dxdy. 



\y\<r 
\x\>(l-X)r 



Moreover, we only increase the last integral if we integrate over all y G 
Thus 



If Tr €2 \r F (x,y)\ 2 dxdy < J (Vr c2 J 7 HF (x, y)7 HF (y, s)dy) dx 
\y\<r b|>(l-A)r 

\x\>(l-X)r 

I Tr c2 [(7 HF ) 2 (x,x)] dx. 
\x\>{l-\)r 

If we now use that (7 HF ) 2 = 7 HF and that p HF (x) = Tr C 2 [y HF (x,x)] we obtain 
the estimate 
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(58) ±JjTr c2 \r F (x,y)\ 2 V(x,y)dxdy < Ex + Ex [7 r HF ] 

+ Z. J p ^ {x)dx . 

\x\>(l-\)r 

If we combine (54), (57) and (58) we obtain 

£HF [7 HF] > £HF [7 HF] + £ A [% H P] _ & [% H F] + ^ ^_ 1 A _ $ HF_^ ^ 



(0) 



7T 



(A(l - A)r)" 



|x|>(l-A)r 



Since < 7 ( q F < / and the density of 7™ is supported within the set 



(0) 



{x : (1- A)r < |x| < (1 - A) _1 r} 
we have from the Lieb-Thirring inequality (22) that 



Tr 



7, 



(") 



> -2Li 



/ 



(l-A)r 



(X) 



5/2 



(l-A)r<|a;|<(l-A)- 1 r 



The factor of 2 above is due to the spin degrees of freedom. We have thus 
proved the lower bound in (53). □ 

As a consequence of this theorem and the Lieb-Thirring inequality (21) 
we get the following bound. 

Corollary 6.3 (L 5 / 3 bound on p™). Let K\ denote the constant in the 
LT inequality (21) and eo, as in (9), denote the TF energy of a neutral atom 
with unit nuclear charge and physical parameter values. Then 

7/3 



(59) J P? ¥ (yfl*dy < 2K^K + |(37r 2 ) 2 / 3 ifr 2 e 
where TZ was given in (51). 



r sup <P r (X) 

\x\=r 



+ 



Proof. Since <& HF is harmonic on the set {|x| > r} and tends to zero at 
infinity we get for all \y\ > r that <E> r (y) < |y|~ 1 r sup^^,, $ r (x). Hence 



r sup 9 r (x) 

\x\=r 



£ A [ir\ > K x J dy 

+ \jjpf F (v)\v-vr 1 pr{y')dydy' 



J \y\- 1 pf Y (y)dy 



From standard atomic Thomas-Fermi theory it follows that the right-hand side 
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is bounded below by the energy of a neutral Thomas-Fermi atom with nuclear 



charge 
A simp 



rsup N=r $ r (x) 
e scaling argument shows that this is 



and with the constant K\ in front of the first term. 



r sup <I> r (x) 

\x\=r 



7/3 



e . 



By repeating this argument with only a fraction of the term / (P J HF ) 5 ^ 3 we 
conclude that for all < t < 1 



(1 - t)K x J Pr(y) 5/3 dy < £ A [7 r HF ] + U^ 2 ) 2/S ( tK ^~ le o 



1 7/3 



r sup <E> r (x) 

\x\=r 



Since £ A [l^ F ] < 71 (by choosing the trial 7 = in (50) ) we get (59) if we 
choose t = 1/2. □ 

We still need to show how we can control the exchange term Ex [l™] ■ This 
is done using a standard inequality of Lieb [8] (or in an improved version by 
Lieb and Oxford [14]). They proved the inequality for general wave functions, 
but we need it here only for Hartree-Fock Slater determinants or more precisely 
for density matrices. For completeness we shall give a proof (with a worse 
constant) in the simple case we need here. 

Theorem 6.4 (Exchange inequality). For any trace class operator 7 with 
< 7 < / we have the estimate 



£*[y] < 1-68 J pf 3 . 



Proof. We shall here present a simple proof that the inequality holds with 
1.68 replaced by 248.3. To get the much better constant one needs the more 
detailed analysis in [14]. We use the representation 

roo 

|x| _1 = tt' 1 / X r * X r (x)r~ 5 dr, 
Jo 

where X r again denotes the characteristic function of the ball of radius r cen- 
tered at the origin. Thus we may write the exchange energy as 

roo p 

£ X [ 7 ] = (27T)- 1 / / TT^Xr^Xr^dzdr, 

Jo Jm. 3 

where X r ^ z is the multiplication operator X r ^ z f(x) = X r (x — z)f(x). 

We now use the two simple estimates X r ^X r ^ z < X^ z and X riZ jX rjZ < 
Tv[-fX^ z ]L We obtain 
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TlfrXr^X^} = Tr^Xr^Xr^Xr^X^J 1 / 2 } < Tr[ 7 X 2 J = Pl * Xr(z) 

and 

Tr[ 7 A^ 7 X r , z = Trtf^X^iXr^Xr^Xr^ 1 / 2 ] < Tr[ 7 X r %] 2 = (p 7 * Xr(z)f . 

If p* denotes the Hardy-Littlewood maximal function we have Pj*Xr(z) < 
(Airr 3 /3)p*(z). Thus with R(z) = ((47r/3)p*(z))~ ' we can estimate 

£x [ 7 ] < (2.)- £ (/ o ^_ W j r-dr + — W r-*-J * 

= ( 4 -/3) 1/3 / R3 (p;w) 4/3 ^. 

If we finally apply the maximal inequality ||/*||£ < (48p2 p /(7r(p- 1))||/||£, 
for all p > 1 ([28], p. 58) we obtain 



7. Exterior L -estimate 

The aim of this section is to control Jj x | >f . P KF (x)dx, for all r > 0. As 
before p HF is the density of a Hartree-Fock minimizer 7 HF with Tr[ 7 HF ] = N. 
Thus / P HF = N. 

The difficulty in estimating Ji x i >r p HF (x)dx is that this quantity cannot 
be controlled in terms of the energy £ HF ( 7 HF ). More precisely, Jj x | >r P 1 {x)dx 
can be arbitrarily large even when £ HF ( 7 ) is arbitrarily close to the absolute 
minimum. The simple reason is that "adding electrons at infinity" will not 
raise the energy. 

Therefore, in order to control j\ x \ >r P HF (x)dx, we must use the minimizing 
property of 7 HF . 

In contrast, it follows from the Lieb-Thirring inequality that 
I\x\>r P HF (x) 5 ^dx can be controlled in terms of the energy. By Holder's in- 
equality it then also follows that the integral of p HF over any bounded set can 
be controlled by the energy. 

The philosophy here will be, to use the minimizing property of 7 HF , to 
control the integral of p HF over an unbounded set, in terms of the integral over 
a bounded set. 

Our main result in this section is stated in the next lemma. The proof of 
the lemma uses an idea of Lieb [13]. 
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Lemma 7.1 (Exterior L 1 -estimate). For all r > and all < A < 1 the 
density p HF of an HF minimizer 7 HF satisfies the bound 



( P HF (x)dx < 1 + 2X^ + 2 

^IxXl-A)-^ 



sup |x|$ (1 _ A)r (x) 

|x|=(l-A)r 



+ 

/ \ 1/2 

+ [K x r- 1 [ p™{x)dx) , 

\ Jr<|x|<(l-A)- 1 r / 



w/iere := (^ + (1 — A) x ) (^x) 2 - -Here < I ) | ! 1 F _ A - )r is £/ie screened nuclear 
potential introduced in Definition 1.1. 

Proof. Since 7 HF is a minimizer we know that it satisfies the Hartree-Fock 
equations. For example, according to Theorem 3.11, 7 HF is a projection onto a 
space spanned by functions u±, . . . , ujv £ L 2 (R 3 ; C 2 ) satisfying H^hfUi = SiUi, 
were E{ < 0. 

Let S G C^M 3 ) have compact support away from x = 0, be real and 
satisfy H(x) 2 < 1. Then 



> J>/ M x)| 2 |x|~(x) 2 dx = £ /^(xH^^/^^x. 

i=l i=l 

From the definition (30) of the mean field operator H^hf we obtain 

N 

(60) > / V(u;(x)*|x|~(x) 2 ) • Vu;(x)dx-Z / p HF ~ 2 

|2/I^(y) 2 



1=1 

+ 



p HF (x)p HF (2/) - Tr c2 |7 HF (x,y)r 



F ~ 2/1 



-dx dy. 



We consider separately the different terms above. By the IMS formula (19) 
(61) Re± J V (ui(x)*|x|~(x) 2 ) • Vu;(x)dx 

= i y | V (-(x)^ 1 / 2 ^)) | 2 - IM^-Or) 2 k(x)| 2 dx 
-i | (IV (H(x) 2 ) + |x| (VH(x)) 2 ) Mx)| 2 dx 

> ~\ j (l V ( S ^) 2 ) + 1*1 (VS(x)) 2 ) h(x)| 2 dx, 

where we have used Hardy's inequality / |V/(x)| 2 dx >\f |x| _2 |/(x)| 2 dx. 
For the Coulomb terms we estimate it using that 

N 

p HF (xK F (y)-Tr c2 [|7 HF (x,y)| 2 ] = \ £ \\ Ul {x) ® u 3 {y) - u 3 {x) ® Ul {y)\\l W 
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is nonnegative. Hence 



P HF (x)p HF (y)-Tr c2 |7 HF (x,y)| 2 



= |/[p HF (xK F (y)-Tr c2 [| 7 HF (x, 



\y\-{yf 



\x - y\ 



dx dy 



\y\(l-E(x) 2 )E{yf 



dx dy 



+ 



p HF (x)p™(y)-Tr C 2 |7 HF (x,y)| 2 



\x-y\ 
\x\ + |y|„ 2 2 



\x - y\ 



-E(x) E(y) dxdy, 



where we expressed the last term symmetrically in x and y. If we now use the 
triangle inequality and the fact /Tr C 2 |7 HF (x, y)| 2 E(y) 2 dy < p HF (x), which 
follows from E(x) 2 < 1 and (7 HF ) 2 = 7 HF , we arrive at 



(62) 



/ / \p HF (x)p UF (y) -Tr c2 [|7 HF (x,y)| 2 ll f^-dxdy 

-II [p HF ( x )p llF (y)- T ^[\^ F (x, y )\ 



\y\(l-E(xy)E( y ) 
\x - y\ 



■dx dy 



+ 1 ( / pHK H 2 j _ 1 j p HP S 2_ 



Inserting the inequalities (61) and (62) into (60) gives 



(63) 



> -\ j (±V (e(x) 2 ) + \x\ (VE(x)) 2 ) p nF (x)dx -Z J p HF E 2 

|y| (l - E(z) 2 ) E(y) 2 



+ 



+ 



/> HF (*K F (y)-Tr c2 |7 HF (x,y)| 



\x-y\ 



dx dy 



p HF E 2 



1 / pHF"2 



By an approximation argument it is clear that we can use (63) for any 
real function E for which E 2 < 1 and the function (\7 (E(x) 2 ) + \x\ (VS(s)) 2 j 
is bounded. In particular we can choose S identically equal to 1, which will 
recover Lieb's result from [13], i.e., / p HF <2Z + 1. 

We shall now choose E := 6 r , where r is the localization function given 
in Definition 6.1. Then 



(64) IV (E(x) 2 ) + \x\ (VE(x)) 2 < JL + (1 - A)-V^ = K x r^ 
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and 



J f [p UF (x) P ™(y) -Tr c2 [| 7 HF (x,y)| 2 ]] M (l ^ {v? dx dy 

> 

-dx dy. 



p HF (xK F (y)-Tr c2 |7 HF (x,y)| 2 



y\ 

\y\O r (y) 2 

\x-y\ 



|x|<(l-A)r 

If we now use that \x\ < (1 — X)r and y G supp 8 r imply that — < A" 1 
we obtain 

|y|(l-~(z) 2 )~(y) 2 



(65) J f [p HF (x)p HF (y) - Tr c2 [| 7 HF (*,y)| 



\x-y\ 



-dx dy 



> J J p^( x )p^{y)\^Ldxdy - A- 1 J e r (y) 2 p^(y)dy, 
\x\<(l-X)r 

where we have also used that jTr^i |7 HF (ir, y)| 2 dx = p HF (y). 
If we insert (64) and (65) into (63) we arrive at 



> [ P HF (x)dx- fe r (yf\y\^ (y)p^(y)dy 

*r Jr<\x\<(l-X)- l r J ( ' 

+ \(Jp™eif-(\ + ^)j P ™el 



Using now that ^^_ x y{y) tends to zero at infinity and is harmonic for \y\ > 
(1 — A)r, which contains the support of 8 r , we see by a simple comparison 
argument that 



9r(yy\y\^_ x)r (y)<e r ( y y 



sup \x\$ (1 _ x) (x) 

Jx|=(l-A)r 1 > 



■J + 



Thus 

> 



-,HF/)2 



- 1 + 2A~ X + 2 



sup |x|$ (1 _ A)r (x) 

|>|=(l-A)r 



nHF/]2 



-Kxr- 1 f p HF (x)dx. 

Jr<\x\<(l-\)- 1 r 

Finally, in order to arrive at the result of the lemma we simply use that > 
X 2 - BX - C for B, C > implies X < B + y/C. □ 



8. The semiclassical estimates 

In this section we derive the relevant semiclassical estimates. We do not 
attempt to give optimal results. We shall be satisfied with what is needed for 
the application we have in mind. In a certain sense it is misleading to refer 
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to the estimates in this section as semiclassical. Usually, semiclassics refers to 
the limit as Planck's constant Ti tends to zero. One then expands the relevant 
physical quantities like energy and density in powers of fi. In our setting there 
is, however, no semiclassical parameter which could play the role of Planck's 
constant. It is rather that we consider potentials for which the semiclassical 
expressions for the energy and density are approximately valid. We must then 
estimate the errors directly in terms of certain norms of the potential. The 
estimates are semiclassical in the sense that if one introduces a semiclassical 
parameter then the errors are of smaller order than the leading semiclassical 
expression. 

We are interested in a semiclassical approximation to the sum of the neg- 
ative eigenvalues of a Schrodinger operator 

h := -£A- V, 

on R 3 . We shall in this section always assume that the potential V : R 3 — > R 
is locally in L 1 and that its positive part satisfies that [V]+ G L 5 / 2 (R 3 ). This 
ensures (by the Lieb-Thirring inequality or even by Sobolev's inequality) that 
h is bounded below and can be defined as a Friedrichs' extension from the 
domain Cg°(R 3 ). 3 

The semiclassical approximation to the sum of the negative eigenvalues of 
h is given by 

(66) (2vr)- 3 J J \p 2 - V(x) dpdx = -2 3 / 2 (15vr 2 )- i y [V»]+ /2 dx. 

\p 2 ~V{x)<0 

Moreover, the semiclassical approximation to the density, i.e., the sum of the 
absolute square of the eigenfunctions corresponding to the negative eigenvalues, 
is 

(67) (2vr)- 3 f ldp = 2 3 / 2 (6tt 2 )- 1 [V(x)]f . 

J±p 2 -V{x)<0 



Definition 8.1. For s > 0, let g : R 3 — ► R be the ground state of the 
Dirichlet-Laplacian for a ball of radius s, i.e., the function with g(x) = if 
\x\ > s and g(x) = (27rs) _1 / 2 |a:| _1 sin(7r|x|/s) if \x\ < s. Then 

(68) < g < 1, Jg 2 = 1 J \Vg\ 2 = (tt/s) 2 . 



3 Notc that wc are not including spin. The operator h is acting in the space L 2 (R ) and not in 
the space L 2 (M 3 ;C 2 ). 
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Lemma 8.2 (Semiclassical approximation). We assume about the poten- 
tial that [V]+, [V - V * g 2 ]+ G L 5 / 2 (R 3 ) with g as in Definition 8.1 above. Let 
e (!) < e ( 2 ) < . . . < 0, denote the negative eigenvalues of h = — ^A — V as an 
operator on L 2 (M. 3 ). Then for all < 5 < 1, all integers N > 0, and all s > 
we have 

(69) 



N 

£ e (j) > — 2 3 / 2 (157r 2 ) _1 (l — <5)~ 3 / 2 / [V] 



5/2 
+ 



■±ttV 2 ^--M~ 3/2 



F - F * 5 " 



5/2 
5/2 



where L\ is the constant in the Lieb-Thirring estimate (22). If there are fewer 
than N negative eigenvalues the sum on the left refers simply to the sum over 
all the negative eigenvalues. 

On the other hand, if also [V]+ £ L 3 / 2 (R 3 ), me can, for all s > 0, /md a 
density matrix 7 mi/i p 7 (x) = 2 3 / 2 (67r 2 )~ 1 [V]+ 2 * g 2 {x) such that 



(70) 



Tr[-lA 7 ] = 2 1 / 2 (5vr 2 )- 1 | [F]f + \^s~ 2 f 2 3 / 2 (6^ 2 )- 1 [V] 



3/2 
+ 



Remark 8.3. We are not proving that the true density of the projection 
onto the negative eigenvalues of h is approximated by the semiclassical expres- 
sion (67). We only claim that there is a 'good' trial density matrix. In the 
context where we shall use the semiclassics we shall infer the approximation of 
the true density by other means. 

Proof of Lemma 8.2. We prove the result using the method of coherent 
states (see Thirring [29] and Lieb [10]). 

For u,p £ M 3 let n up be the one-dimensional projection in L 2 (M 3 ), pro- 
jecting onto the space spanned by the function f u ,p(x) ■= exp(ipx)g(x — u). 
We then have the coherent states identities 

(71) Tr[n UjP ] = 1, for all p, u 

(72) (2vr)" 3 J J U UtP dpdu = /, on L 2 (R 3 ). 
We also have the identity 

(73) Tr[-lAn u , p ] = \p 2 + \J \Vg\ 2 = \p 2 + n 2 /(2s 2 ), 
and for all density matrices 7 

(74) Tr[-iA 7 ] = (27r)- 3 ||ip 2 Tr[n UiP 7]dpdn-7r 2 /(2 S 2 )Tr[ 7 ] 

(75) Tr[(y*<7 2 ) 7 ] = (2vrr 3 J J V (u)Ti[U u ^]dp du . 
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Proof of the lower bound (69). Let f±, fa, . . . G L 2 (M 3 ) denote the normal- 
ized eigenfunctions of h corresponding to negative eigenvalues. It is clear that 
we may without loss of generality assume that there are N negative eigenvalues. 

We decompose the operator h as 

h = - 6) A - V * g 2 + [-\8A - {V - V * g 2 )]. 
We may then write £)JLi = EjLi(/j, hfj) =A + B, where 

N N 
j=l 3=1 

Hence from (74) and (75) we obtain 

r r N 
A = (27T)- 3 / / (j(l - 5)p 2 - V(uj) ^(fj^u, P fj)d P du - NTt 2 {2s 2 )-\ 

J J j= i 

As a consequence of (71) we have that < J2jLi(fj^u,pfj) < 1- 
It is therefore clear that 

A > {2tt)- 3 JJ ({l-5)\p 2 -V{u))dpdu-T: 2 {2s 2 )- 1 N 

(l-5)ip 2 -V(«)<0 

= -2 3 / 2 (15vr 2 )- 1 (l - ,5)- 3 / 2 J [V]f - it 2 {2s 2 y l N. 

The estimate (69) follows by applying the Lieb-Thirring estimate (22) to 
conclude that B > -L^^W [V - V * g 2 ] + f^. 

Proof of the existence of 7. We shall prove that 

7 :=(27r)- 3 JJ U U)P dpdu 

ip 2 -V(u)<0 

has the desired properties. From (72) we see that 7 is a density matrix, i.e., 
< 7 < /. The density corresponding to 7 is easily computable 

P 7 (s) =-y(x,x) = (2tt)- 3 JJ n u , p (x,x)dpdu = 2 3 / 2 (6ir 2 y 1 [V} 3 / 2 *g 2 (x) 

ip 2 -V(u)<0 

From (73) we immediately obtain (70). □ 

Although we shall use the semiclassical approximation in the form given in 
the lemma we shall for completeness state a less technical semiclassical result 
which follows very easily from the lemma. 
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Theorem 8.4 (Semiclassical approximation). Assume that < V G 
L 5 / 2 (K 3 ) ni 3 / 2 (l 3 ) and |W| G L 5 / 2 (i? 3 ). Let e« < e( 2 ) < . . . < 0, denote 
i/te negative eigenvalues of h = — - V os an operator on L 2 (K 3 ). T/ten we 
have 

5 

(76) ]Te (i) > -2 3 / 2 (15vr 2 )- i y V» 5/2 dx |l + A L ||F||7* ||W ||| ||F|| | J 3 
and 

(77) Ve^ < -2 3 / 2 (15vr 2 )- 1 / V» 5 / 2 dx + 2" 1 / 2 vr- 4 / 3 1 1 1 1 5 1| W||| 

1 /3 

where A L := §2- 9 / 10 (15tt 2 ) 3 / 5 (^) l£ /3 l} /15 . #ere L and Li are the 
constants in the CLR and Lieb- Thirring estimates (23) and (22) respectively. 



Proof. We may estimate 
- V *g 2 {u) 



\VV(u-ty)\\y\g(y) 2 dtdy 

= f \VV(u-y)\\y\ [ l t- A g{y/t) 2 dtdy 
Jm? Jo 

< (47T)- 1 / \VV(u-y)\\yr 2 dy, 

J\v\<s 



where we have used that t- 4 g(y/t) 2 dt < |y|~ 3 j* °° t 2 g{t) 2 dt = (47r) -1 |j/| -3 
(identifying g with a function on K+). Hence 

(78) ||^-y* ff 2 || 5/2 < (47r)- 1 ||VF|| 5/2 / |yr 2 ^ = s||vi/|| 5/2 . 

•%l<* 

For any density matrix 7, Trf/ry] is an upper bound to the sum of the negative 
eigenvalues of h. From (70) we find for the density matrix constructed in 
Lemma 8.2 that 

i5/2 



Tr[/i 7 ] = -2 3 / 2 (15vr 2 )- i y [V] 

• 2 3 / 2 (6vr 2 r 1 J [V{u)fl 2 [V(u) - V * g 2 (u) + ±vr 



+ 



2 s' 2 



du. 



The bound (77) follows from applying Holder's inequality, (78), and optimizing 



m s. 



By the CLR bound (23) we know that h has only finitely many negative 
eigenvalues and that their number N is bounded by N < Lq f [V] 3 / 2 . From 
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(69) and (78) we therefore obtain 
£ e (i) > -2 3 / 2 (15tt 2 )- 1 (1- 5)- 3 / 2 Jv 5 / 2 

3 

-L VT 2 (2 S 2 )- 1 J ^3/2 _ Lir 3/2 s 5/2 J | W |5/2 

= -2 3 / 2 (15tt 2 )- 1 (1-<5)- 3/2 J V 5 / 2 

- 1 (|) 5/9 4 /9 ^ /9 - 10/9 ^ 2/3 (/ ^ 3/2 ) 5/9 (/ iwi 5 / 2 ) 479 , 

where we have optimized in the parameter s. 

We now optimize in the parameter S. Define 5' by (1— 5)~ 3 / 2 = (1— 5')~ 2 ^. 
(Note that < 5 < 1 if and only if < 8' < 1.) Then <T 2 / 3 < (4y/9)~ 2 / 3 . 
Thus 

ax 5/9 / , x 4/9 

F 3 / 2 ] ( / |W| 5 / 2 ] , 

where A 1 := (^-f L f 'if \ Using that 

min (1 — 

,5')-2/3 + < y'-2/3 6 j = o[1 + ( 6 / a )3/5]5/3 

we arrive at (76). □ 

We shall need the semiclassical estimates also for the operator h restricted 
to functions on the set {x : \x\ > r} satisfying Dirichlet boundary conditions. 

Lemma 8.5 (Dirichlet boundary conditions). Let the assumptions be as 
in the beginning of Lemma 8.2. For r > let h r denote the restriction of 
the operator h = — |A — V to functions on the set {x : \x\ > r} satisfying 
Dirichlet boundary conditions. Denote by < e^ < . . . < the negative 
eigenvalues of h and by ef 1 ^ < e? < . . . < the negative eigenvalues of 
h r . Then Y^j e ^ — J2j e r^- Moreover, if 7 is a density matrix on L 2 (IR 3 ) we 
may, for all < A < 1, find a density matrix 7 such that p^ is supported in 
{x : \x\ > r} and p^ < p 7 and 

Tr[h r j] < Tr[H + L 1 f [V\% 2 + i(vr/(2Ar)) 2 f P T 

J|x|<(l-A)- 1 r ^Ixl^Cl-A)-^ 



Proof. That the Dirichlet eigenvalues are upper bounds to the eigenvalues 
on M 3 , is a well known simple consequence of the variational principle. 
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Let 9 r be the localization function from Definition 6.1. We shall choose 
7 = e r <y6 r . Then by the IMS formula (20) we have 

Tr[/i 7 ] = Tr[h r 9 r jd r ] + Tr[/i(l - # 2 ) 1/2 7(1 - #r) 1/2 ] 

-lTr[((V^) 2 + (V(l-^ 2 ) 1/2 ) 2 )7]- 

By the Lieb-Thirring inequality (22) we have Tr[/i(l - 6»^) 1 / 2 -y(l - 9 2 ) 1/2 ] > 
—L\ Jj x |<(i_^)-i r [V]+ 2 ■ Thus the lemma follows from the bound on the gra- 
dient of 6 r and (1 - 0j!)i/2 

given in Definition 6.1. □ 



9. The Coulomb norm estimates 

In this section we introduce and study the Coulomb norm. 
Definition 9.1. For /, <? G L 6 / 5 (R 3 ) we define the Coulomb inner product 
(79) D(f,g) :=\JJ f {x)\x - y^gl^dx dy 

and the corresponding Coulomb norm, 

\\g\\ c :=D{g,g) l l 2 . 



By the Hardy-Littlewood-Sobolev estimate we have 

(80) lbllc<vr- 1 /6 2 V6 3 -i/2 ||5||6/5 . 

In this sharp form the inequality was proved by Lieb [12]. Using the Fourier 
transform we may write 

(81) D(f,g) = (2tt) J f(p)W)\p\- 2 dp, 

from which it follows that the Coulomb norm really is a norm on L 6 / 5 (M 3 ). 

The following estimate was first used in the context of atomic problems 
by Fefferman and Seco [5]. 

Lemma 9.2 (Coulomb norm estimate). /// G L 6 (M 3 ) with |V/| G I? 
and g G L 6 / 5 (M 3 ) then 



fg 



< (2vr)- 1 /2|| V /|| 2 || 5 || c . 



Proof. Using Plancherel's identity and the representation (81) we have 

/ fg\ = |/ fl\ < ||N/(p)|| 2 ||br^(p)||2 = (2^- 1 / 2 ||V/|| 2 ||< 7 || c . □ 
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We shall next give some simple but very useful consequences of this esti- 
mate. 

Corollary 9.3. Consider f G L 5 / 3 (R 3 ) n L 6/5 (M 3 ). For all x£R 3 and 
s > we have 

(82) /* \x\- 1 < (25^ S /16) 1 /5 (vrs )i/5|| [/]+ || L5/3(B{ ^ )) + (2/ s ) 1 /2|| / || c . 

For all x G IR 3 and all k > denote by A(\x\,k) the annulus 
A(\x\,K) = {y:(l-2K)\x\ < \y\ < \x\} 

we then have 

(83) f f(y)\x-y\' l dy < 2 7 /V/ 5 (kM) 1 /5 || [/]+ || 
J\y\<\x\ 

+ 2 3 / 2 «r 1 |*r 1/2 ll/Hc. 



Remark 9.4. Note that we do not restrict to k < 1/2. We do this to 
avoid having to check this condition in the applications of the corollary. 



Proof of Corollary 9.3. Consider the function £ s : 



defined by 



if \z\ < s 
if \z\ > s. 



It satisfies ||V£ S ||2 = (^ir/s) 1 / 2 . Hence from Lemma 9.2 we obtain 

/*N _1 < / [f(y)] + (\x-y\' 1 -s- 1 )dy+ [ f(y)Ux-y)dy 
J\y-x\<s v 7 JR 3 

< (25 7 r 4 S /16) 1 / 5 ||[/] + || L 5/3 (B(a;iS)) + (2/ S ) 1 / 2 ||/||c, 



5n^ 



where we have used that Jj^^idz/I 1 — l) 5 ^ 2 dy - , . 

In order to prove the second half of the corollary we introduce the function 
E X)K : R 3 R 3 given by 



^xA z ) '■= ' 



1, if \z\ < (1 - 2k)|x| 

1 - OIk^OI - |x|(l - 2k)), if (1 - 2k)\x\ < \z\ < (1 - k)\x\ 
0, if (1 - k)\x\ < \z\. 



Then |S XjK (2;)| < 1 and we can estimate 

/ Mlx-Vl^dy < I [fiyMx-y^dy 

J\y\<\x\ JA(\x\,k) 

+ / f{v)\x - y\~ lr Ex,K{y)dy 

JR 3 
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(\ -Vt> 
Lj x - yr5/2dy ) 

+ (2vr)-V2 ||/|| c y | V , (\x - vr^iv)) f dy 



1/2 



It remains to estimate the two integrals above. For the first integral we find 
for k < 1/2 

/ \x-y\~ 5 / 2 dy = 2tt\x\ 1 I 2 ( l (1 - 2ru + r 2 y 5 / 4 r 2 dudr 

JA(\x\,k) Jr=\-lKJu=-\ 

= 2^|x| 1 / 2 oi(k), 

where 

ai(«) := 4(2k) 1 / 2 - (4/3)(2k) 3 / 2 + (4/3)>/2[l - (1 - k) 1/2 (1 + 2k)] < 4(2k) 1 / 2 . 

The last inequality follows from a straightforward careful analysis of ai(re). 
For the second integral we get 



/ V tf (|z - yr^fo)) "dy) < (/ |x-y|- 4 ^) 

'J \J|x-j/|>k|x| / 



\ 1/2 



with 

oli(k) := 1 + (1 — ft) In 
where we have used that 



+ (k\x\) 1 [ / \x-y\ 2 dy) 

\JA(\x\,k) J 
= (47T) 1 / 2 ( K | x |)-V2 (1+a2(K) l/2 ) 

1 — ft\ 1 

< 1 — m k < k , 



\ 1/2 



/ \x— y\ 2 dy = 2ir\x\ \ / (1 — 2rn+r 2 ) 1 r 2 du dr = 47tk|x|q;2(^)- 

JA(\x\,k) J r =l-2n Ju=-1 

Using (4tt) 1 / 2 k- 1 / 2 (1 + q 2 (k) 1/2 ) < 4tt 1 / 2 k- 1 we get (83). 

The estimate holds also for k > 1/2 since the last term in (83) can be 
ignored in this case. □ 



10. Main estimate 

We now restrict attention to the case N > Z. Throughout the remaining 
part of this paper p TF , ip T¥ , <3?^ F , p HF , </? HF and $™ always refer to the problems 
with particle number N. In fact, since N > Z the TF functions correspond 
to the neutral atom, i.e., /U TF = 0. We shall suppress the dependence on N 
everywhere since it is held fixed throughout the discussion. 
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From now on we shall no longer explicitly compute the constants involved 
in the estimates. We shall use the notation (Const.) to refer to any universal 
(in principle explicitly computable) positive constant. Thus (Const.) does not 
mean the same constant in all equations or inequalities. Even within the 
same equation we shall use the notation (Const.) to refer to possibly different 
universal constants. Universal constants of particular importance will be given 
separate names. We begin by stating the main result of this section. 

Theorem 10.1 (Main estimate). Assume Z > 1 and N > Z. There 
exist universal constants < e < 4 and Cm, C$ > such that for all i£l 3 
we have 



(84) *™(x) - *™(x) < C*|x|- 4+£ + C 



'M- 



We shall prove Theorem 10.1 by an iterative procedure. The first step is 
to control "small" x. 

Lemma 10.2 (Control of the region close to the nucleus). Assume Z > 1 
and N > Z. For all (3 > and all \x\ < f3Z~ 1 / 3 we have 



iHF/ \ iTP/ i 



< A$, 



where e\ = 1/66 and A$ > is a universal constant. 



12-£ll |-4+£l 



Lemma 10.3 (Iterative step). Assume N > Z. For all 5,e',a > with 
5 < 5q, where 5o is some universal constant, there exists constants £2,C$ > 
depending only on 5 and a constant D = D(e', a) > depending only on e', a 
with the following property. For all R$ < D satisfying that PqZ^ 1 ^ < Rq +S 

(9tt) 2 / 3 



(where (5q 



44 



as in Theorem 5.4) and that 



(85) 



< a\x 



-4+s' 



holds for all \x\ < Rq, there exists R' > Rq such that 



(86) 



< CL\x\- 4+£2 



for all x with Rq < \x\ < R' Q . 



Lemmas 10.2 and 10.3 will allow us to control small and intermediate \x\ 
to control large |x| we shall need the following two lemmas. 
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Lemma 10.4 (Bound on J(p HF ) 5//3 )- Assume N > Z. Given < e',a, 
there is a D > such that if (85) holds for all \x\ < D then we have 

(87) / P* F (y) 5/3 dy < (Const.)|x|- 7 , 

J\y\>\x\ 

for all \x\ < D. 



Lemma 10.5 (Bound on f p HF ). Assume that (85) holds for all \x\ < R 
for some R > and some e' , a > 0. Then for < r < R we have 



\y\<r 



(p HF (y)-p TF (y))dy 



< o-r~ 3+£ ' 



and 



[ P HF (y)dy < (Const.)(l + or £ ') (r~ 3 + 1 
J\y\>r V 



We shall prove Lemma 10.2 in Section 11 and Lemmas 10.3 and 10.4 in 
Section 12. We end this section with the proofs of Lemma 10.5 and the main 
estimate Theorem 10.1. 

Proof of Lemma 10.5. First note that for < r < R we have 
/ ( P TF (y)-p HF (y))d y =(4n)~ 1 r f f ( P ^ F {y)-p^ F (y))\ru-y\~ 1 dydu. 

J\y\<r J J 



ui£S 2 \y\<r 



Thus we have 



/ (p TF (y)-p ilF (y))dy=(A7rr 1 r [ ^ F (ru) - ^ F (ru)du. 

J\y\<r J 

uies 2 

Together with (85) this gives the first estimate above. Moreover, we also have 
that 



(P TF (y)-P HF (y))dy 



< sup 

\y\=r 



C F (y) - C F (y) 



lr/2<\y\<r 

+ sup &(y)-&(y) 

\y\=r/2 
< (Const>r" 3+£ '. 

The TF equation (6), and the Sommerfeld estimate in Theorem 5.2 give 

/ P TF (y)dy < (Const. )r~ 3 
J\y\>r/2 

andhenCe / 

Jrl 



r/2<\y\<r 



P HF (y)dy < (Const.)(l + ar £ )r" 3 . 
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From (85), the exterior L 1 -estimate Lemma 7.1 (used with A = 1/2 and r 
replaced by r/2), and Lemma 5.5 (recall that now ^ TF = 0) we immediately 
conclude the estimate on f\ y \ >r P R¥ (y)dy. 



□ 



We finally show how to use Lemmas 10.2-10.5 to prove the main estimate 
Theorem 10.1. 

Proof of Theorem 10.1. We first show that we may choose 5 > small 
enough such that if we choose R 1+5 = PqZ- 1 / 3 we have for all \x\ < R that 

(88) *™{x)-*Z(x\ <C'i\x\-^ 



for a universal constant C$ > and with E\ given in Lemma 10.2. 

To see this let (3 > be such that (/?Z -1 / 3 ) 1+ * = A)^~ 1/3 , i-e., 
/3 Z 5 / 3 . We then see from Lemma 10.2 that for all \x\ < j3Z~ 1 / 3 we have 



i+s _ 



*jj(s) - *j» 



° x 6 \x\ ' 



< ~Z~ 

= A#(P Z s / 3 )(t5^-2&v)z- ! &-\x\- 4+ %. 

Since R 1+s = PqZ^ 1 / 3 , i.e. R = f3Z~ 1 / 3 , we see that if 5 is small enough and 
C$ is chosen appropriately then (88) holds for all \x\ < R. 

We now assume that 5 is also small enough that we may apply Lemma 10.3. 
This gives us constants £2)C$ > and for all a, e' > a D > with the 
properties stated in Lemmas 10.3 and 10.4. We may without loss of generality 
assume that D < 1. Now choose a = max{ C'$,C§} and e' = min{ei/2, £2}- 
Note that a, e, and D are now universal constants. We shall prove that for all 
\x\ < D we have 



(89) 



,HF/ n 



< o~\x\ 



-4+e' 



Since we are assuming that D < 1 it is sufficient to prove (89) with e' replaced 
by ei/2 or £2- We have to prove that D belongs to the set 

M = {0 < R < 1 : Inequality (89) holds for all \x\ < R}. 

If this were not true we would have D > Rq := sup.M. In order to reach a 
contradiction we therefore assume this and hence in particular that i?o < 1. 
From (88) it follows that either R > 1 or R £ M. If R > 1 then R = supM 
= 1 which contradicts our assumption. On the other hand if R £ M then 
i?Q +<5 > R 1+s = PqZ~ 1 / 3 . It is then an immediate consequence of Lemma 10.3 
that there exists R' £ M. with R' > Rq and this is of course also a contradic- 
tion. This establishes an inequality of the form (84) for all \x\ < D. 
We shall now prove (84) for |x| > D. We write 



iHF / \ iTF/ \ 



< 



+ 



D<\ v \<\x\ \x-y\ 



dy 
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We shall estimate the last term using Lemma 10.4 and the similar bound 

P TF (yf /3 dy < (Const.)|xr 7 , 



/ 

J\v 



>\x\ 



which holds for all x by the Sommerfeld estimate Theorem 5.2 and the TF 
equation (6). Hence using Holder's inequality we have 



(p TF (y)-P llF (y)) 



D<\y\<\x\ 



\x - y\ 



dy 



x — y\ 5 ^ 2 dy 
\<b 



< (Const.)!)- 21 / 5 ( [ 

\J\x-y 

+D- 1 I (P TF (y) + P HF (y))dy. 

J\v\>D 



2/5 



By Lemma 10.5 and the bound S\ y \>D P TF (v) dy < (Const. )D _3 , which is again 
a consequence of the Sommerfeld estimate Theorem 5.2 and the TF equation 
(6), we see that this last expression is bounded by a universal constant. 

Since 3> D (x) — <& D (x) is harmonic for \x\ > D and tends to zero at infinity 
we have for all \x\ > D that 



^ F (x) - <b TF (x) < sup ^ F (z) - d> D (,: 
\z\=D 



■TP , 



which is also bounded by a universal constant. Thus (84) holds for all x. □ 



11. Control of the region close to the nucleus: proof of Lemma 10.2 

In order to prove Lemma 10.2 we need some basic estimates. 

Lemma 11.1 (Global L 5//3 and Coulomb norm estimates). For all N and 
Z we have the bound 

(90) / P HF (yf /3 dy < (Const.)Z 7 / 3 . 

Moreover, if Z > 1 
(91) 

with e 3 := 2/77. 



| /9 hf_ p tf||2 < (Const.)^ 1 ^ 3 )/ 3 , 



Proof. Although we shall only use this result for iV > Z the proof is 
almost as easy without this restriction, so we treat the more general case here. 

We first estimate the L 5 / 3 norm of p HF . It is easy to see that £ HF (7 HF ) < 0. 
Thus since P(7 HF ) — <£*(7 HF ) > we have from the Lieb-Thirring inequality 
(21) that 

0>£ HF (7 HF ) > TV[(-iA-Z|y|- 1 )7 HF " 

> J (K lP « F { y fl* - Z\y\- l p« F (y))dy. 
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If we use that / p HF = N and the inequality ab < |a 5//3 + §6 5//2 we get for all 
5 > and all r > 0, 

> h Kl - l5 5 / 3 )p™(yf 3 dy - | [ {5- l Z\y\- l f 2 dy - NZr' 1 . 
J J\y\<r 



Choosing (5 5 / 3 = 5K 1 /6 and optimizing in r gives / (p HF ) 5 / 3 < (Const. )N 1 / 3 Z 2 . 
If we use that, since there exists an HF minimizer with particle number N, we 
must have Lieb's bound N < 2Z + 1 (see Theorem 3.5) we arrive at (90). 

We turn to the proof of (91). We rewrite the Hartree-Fock functional (26) 

as 

(92) £ HF ( 7 ) = Tr [(-IA - ^ F ) 7 ] + ||p TP - pj 2 c - D(p™,p™) - &r( 7 ), 
where we have used the definition tp TF (y) = Z\y\~ x — p TF * and 



Tr 



(P TF * lyl" 1 )^ = 2Z)(p TP ,p 7 ) = 2?( 7 ) + D(p TF ,p TF ) 



\n TF n II 2 
\P — P-y\\c- 



From the semiclassical estimate (69) and the fact that, when 7 is a density 
matrix with Tr[ 7 ] = N and h is a self-adjoint operator, then Tr[/i 7 ] is an upper 
bound on the sum of the N lowest eigenvalues of h, we find 



Tr 



> -2 5 / 2 (15vr 2 )- 1 (l - 5)- 3 / 2 J [<p TF - p TF ] 
- \it 2 s- 2 N - 2LkT 3 / 2 



5/2 
+ 



TF TF 2 

cp -ip *g 



+ 



for all < 5 < 1 and all s > 0. Recall that the function g was given in 
Definition 8.1. Here we have used the semiclassical estimate for the space 
L 2 (M 3 ;C 2 ). The estimate above therefore has an extra factor of 2 in the first 
and the last term compared to (69). Thus 

5 HF (7 HF ) > -2 5 / 2 (15^ 2 )^ 1 (l - <5)- 3 / 2 J [ip TF - p TF tl 2 - p TF N - D(p TF , p TF ) 



+ I|P TF -P HF HC 



I - Itx 2 s- 2 N - 2Li<T 3/2 



TF TF 2 

cp -ip *g 



+ 



5/2 
5/2 



^(7 HF ). 



Since \y\ 1 — g 2 * \y\ 1 > (because the function \y\ 1 is superharmonic) 
we have 



r 



TF 2 

ip *g 



+ 



^ < Z^ 2 \y\~ x - g 2 * [yr 1 ^ < SttZ^s 1 / 2 , 

5/2 5 / 2 



where we have used that |y| 1 — g 2 * \y\ 1 is nonnegative, bounded by \y\ 
and vanishes for |y| > s. If we insert this above and optimize in s we obtain 
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(93) 

£hf (7 h F) > _2 5 / 2 (15vr 2 )- 1 (l - 5)~ 3 / 2 y [^ TF -/x TF ]+ /2 -/x TF 7V-D(p TF ,p TF ) 
+ ||P TF - P uf \\ 2 c - (Const.)^ 6 / 5 ^ 1 ^ 2 - £ X (l w ). 

We choose 5 := ^Z~ 2 / 33 (this is not optimal). Then for Z > 1 we have 5 < 1/2 
and thus (1 - 5)~ 2 / 3 < 1 + (2 5 / 2 - 2)5. Hence 

(94) £ HF (7 HP ) > -2 5 / 2 (15vr 2 )- i y [^ TF - M TF ]^ /2 - fi TF N - D(p TF , p TF ) 

+ ||/0 TF - P HF ||c " (Const.)Z 7 / 3 - 2 / 33 - &(7 Hr ), 
where we have again used iV < 2Z + 1 and the fact that by (90) we have 

2 5 / 2 (157r 2 r 1 J [<p™-f*]*l 2 = § [(3tt 2 ) 2 / 3 ^ J (p TF ) 5/3 l < (Const.)Z 7 / 3 . 

On the other hand, since 7 HF minimizes <? HF among all density matrices 7 
with Tr[7] < N, we can find an upper bound to £ HF (7 HF ) by choosing an appro- 
priate trial density matrix. We choose as trial matrix 7 an operator which acts 
identically on the two spin components. On each spin component we choose 
it to be the density matrix constructed in Lemma 8.2 satisfying (70) with 
V = if TF —fi TF . Note that the Thomas- Fermi equation (6) and the properties of 
7 stated in Lemma 8.2 imply that p 7 = 2 5 / 2 (6^ 2 )- 1 [(p TF - fi TF ] 3 / 2 *g 2 = p TF *g 2 
(the extra factor of 2 compared to Lemma 8.2 is of course due to the spin de- 
generacy). Thus Tr[7] = J p TF < N. 

From (26) and (70) we find, since £x(j) > 0, that 

£ HF (7) < 2 3 / 2 (5vr 2 )- 1 J [^ F - f F ]f+\K 2 s- 2 N-j Z\y\~ l p 7 (y)dy+D(p„ Py ). 



Since J J g(x — z) 2 \z — w\ 1 g(y — w) 2 dw dz < \x — y\ 1 we see that D(p 7 , p 7 ) < 
D(p TF , P TF ). Thus from the definition (2) of (p TF we can write 

S hf (7) < 2 3 / 2 (5vr 2 )- 1 J b TF - ^ TF ]f - J b TF (y) - M TF ] p' FF (y)dy - ^ F N 

-D(p TF ,p™) + ^ 2 s- 2 N + J Z(\yr 1 -g 2 *\y\- 1 )p™(y)dy. 

If we use the TF equation (6), the estimate p TF {y) < 2 3 / 2 (37r 2 )- 1 Z 3 / 2 |y|- 3 / 2 
which follows from the TF equation, and again the facts that — g 2 *\y\~ l 
is nonnegative, bounded by and vanishes for \y\ > s, we obtain after 

optimizing in s 

(95) £ HF ( 7 ) < -2 5 / 2 (15vr 2 )- 1 J [ip TF - P TF ] 5 / 2 - fi TF N - D(p TF ,p TF ) 
+ (Const.)iV 1/5 Z 2 . 
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Comparing (94) and (95) and recalling that £ HF (7) > £ HF (7 HF ) we get 

that 

||pTF _ p HF||2 < (C onst.)^ 7/3 - 2/15 + (Const.)^ 7/3 ~ 2/33 +^(7 HF )- 

If we finally use the exchange inequality in Theorem 6.4 and the estimate (90) 
we see that 

£*(7 HF ) < 1-68 (p HF f 3 y 2 [J (P HF )) V2 < (Const.)iV 1/2 Z 7 / 6 . 
Inserting this above and again using N < 2Z + 1 we arrive at (91). □ 

End of proof of Lemma 10.2. We write 

•%|<|x| 

Using the Coulomb norm estimate (83) we find 
*5(*)-*SH ^ 2 7 / 5 vr 2 / 5 ( K |x|) 1 /5 max {||pT P||L5/3(R3)) ||pH li5/3(R3) j 

+2 3 / 2 K-vr 1/2 iiP HF -p TF iic 

Thus from Lemma 11.1 and the fact that / (p TF ) 5 ^ 3 < (Const. )Z 1 I Z (which can 
be seen for instance from the Sommerfeld estimate Theorem 5.2 together with 
the TF equation (6)) we obtain 

$* F (x) - &w(x)\ < (Const.)( K |x|) 1 / 5 Z 7 / 5 + (Const.)K-Vr 1/2 ^ 7(1 ~ £3)/6 
= (ConstOlxl 1 / 12 ^ 7 /^ 1 -^)/ 36 , 



where the last equality above follows from choosing the optimal value for k. 
Hence if \x\ < /?Z -1 / 3 we have 

The lemma follows since 763/12 = 1/66. □ 

12. Proof of the iterative step Lemma 10.3 and of Lemma 10.4 

We begin by fixing some < r such that (85) holds for all \x\ < r. 
We shall proceed as for the region close to the nucleus, but instead of 
directly comparing HF and TF. We shall introduce an intermediate TF theory. 

OTF 

Namely, the TF theory defined as in Definition 4.1 from the functional £ r : = 

TF 

£ Vr with the exterior potential V = V r given by 



v r (y) -.= Xr(y)K F (y) 



0, if \y\ < r 

®7(y), if \y\ > r. 
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Here again Xr = 1 — X r is the characteristic function of the set {x : \x\ > r}. 
Note that this potential is harmonic and continuous on \x\ > r. Let p° TF 
denote the minimizer for the TF functional, £ r (p), under the constraint 
/ P < / P HF Xr • Denote the corresponding TF potential by 

^ OTF (y) :=v r {y)- P r*\y\- 1 

and the corresponding chemical potential by /U° TF . We shall prove below (see 
Lemma 12.4) that if r is chosen appropriately then p° TF = 0. 

Note that according to the Thomas-Fermi equation (33), p° TF has support 
on the set {y : \y\ > r}. Since V r on the support of p° TF is the potential 
coming from the true HF density for \y\ < r we may interpret p° TF as the 
TF approximation for only the outside region, i.e., \y\ > r, of the atom. The 
notation OTF refers to Outside TF. 

Lemma 12.1 (Preliminary bounds on TF and OTF functions). Assume 
that N > Z then for all y 

V TF (y)< 3 4 2^n 2 \y\' 4 and p TF (y) < 3 5 2-\\y\- 6 . 

For all \y\ > (5qZ~ 1 ^ we have 

</? TF (y) > (Const.)|y|- 4 and p TF (y) > (Const.)|y|- 6 . 

Given e',a > 0, and r > such that (85) holds for all \x\ < r and ar £ ' < 1 
then for all \y\ > r we have 

P r OTF (y) < (Const.)r- 6 and ^ F (y) < \V r (y)\ = |< F (y)| < (Const.)r- 4 . 



Proof. The upper bounds on the TF functions follow immediately from 
the Atomic Sommerfeld estimate Theorem 5.2) and the TF equation (6) if we 
recall that p TF = 0. The lower bounds follow from Theorem 5.4. 

Since <E>^ F is harmonic on the set {\y\ > r} and tends to zero at infinity it 
follows from the assumptions on e',a,r that for all \y\ > r we have 

|^ F (y)| < sup |^ F (z)| < (Const.)r- 4 , 
\z\=r 

where in the last inequality we have used the iterative assumption (85) and 
Lemma 5.5 for the case p TF = and the fact that <E> r > ip TF > (see e.g., 
Theorem 5.4). The inequality (f° TF (y) < \V r (y)\ is trivial from the definition 
ofv? r OTF . 

Finally, from the TF equation (33) we conclude that for all \y\ >r 



Pr TF (y) < (Const.)V r (y) 3/2 < (Const.)r- 6 . 



□ 
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Lemma 12.2 (Preliminary comparison of HF and TF). Assume that N > Z . 
Given e',a > 0, and r > such that (85) holds for all \x\ < r then 



(96) J xt (P^ 



p UF ) < ar- 3+£ '. 



Proof. We have 

/ (P TF (y) - P HF (y)) dy = (4vr)~ 1 r / (V">u,) - ^/(ru,)) du 

J\y\<r JS 2 V ' 

where dio denotes the surface measure of the unit sphere S 2 . Thus according 
to (85) we have 



/ (p TF (y)-P UF (y))dy 

J\v\<r 



\y\<r 

Since / p TF < N = f p UF we have 



< ar- 3+£ '. 



\rP TF - XrP UF ) < I (P HF (y) - P TF (y)) dy < ar- 3+s '. □ 
v 7 J\y\<r 

For \x\ > r we may write 

(97) ^ F (x) - <S>J F (x) =Ai(r,x) + A 2 (r, x) + A 3 (r, x) , 
where 

(98) A!(r,x) = <f^ rp (x)-i P TP (x), 

(99) A 2 (r,x) = [ [P° TF {y)-P TF {y)]\x-y\- 1 dy, 

J\y\>\x\ 

(100) A 3 (r,x) = f [p^{y)-p^{y)]\x- y \- 1 dy. 

Jr<\y\<\x\ 

We turn first to estimating Ai and Ai- Thus we need to control the 
difference between the full TF approximation and the TF approximation for 
the outside region. Our strategy is to first prove that (p° TF (x) and (p TF (x) 
are close on the set {\x\ = r}. An application of the Sommerfeld estimates in 
Theorem 4.6 will then give excellent control on the difference ip° TF (x) — ip TF (x) 
for all \x\ > r. Controlling the behavior on the set {\x\ = r} is difficult and we 
begin with a weak estimate on the difference between p° TF and p TF . In fact, 
we first estimate the difference in Coulomb norm. 

Lemma 12.3 (Coulomb norm comparison of TF and OTF). Assume N>Z. 
Given constants e',a > there exists a constant D > depending only on e',a 
such that for all r with I3qZ~ 1 / 3 < r < D for which (85) holds for all \x\ < r 
we have 

(101) ||p r OTF - XrV F ||^ < (Const.)^-^'. 
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Here again Xr denotes the characteristic function of the set {y : \y\ > r}. 
Moreover, 

(102) M ° TF < (Const.)c7 1/2 r- 4+ T. 

Proof. To prove this we make a perturbation analysis of the TF functional. 
We introduce the perturbation potential 

(103) W(x) = ^ F (x)-^ F (x). 
Then for all \x\ > r 

$7(x) = V r (x) - W(x). 

Note that W is harmonic for \x\ > r and tends to zero at infinity. Hence, since 
we assume that the iterative assumption (85) holds for |x| = r, we have 

(104) sup \W{x)\ = sup \W{x)\ < o-r- 4+£ '. 

\x\>r \x\=r 

We claim that there exist two functions W\, W 2 with supp W\ C {x : \x\ < 3r} 
and supp W 2 C {x : \x\ > 2r} such that W(x) = W\(x) + W2(x) and 

(105) sup \W!{x)\ < sup \W{x)\ < ar~ A+£ ' 

\x\>r \x\=r 

(106) [\VW 2 \ 2 < 4yrr sup \W{x)\ 2 < An<j 2 r- 7+2£ ' . 

J \x\=r 

In order to prove this we let 



F(x) 



if \x\ < 2r 

(|x| — 2r)r _1 , if 2r < \x\ < 3r 

1 if \x\ > 3r. 



Set Wi{x) = (l-F(x))W(x) and W 2 (x) = F(x)W{x). The first estimate (105) 
follows immediately from (104). By a simple integration by parts, similar to the 
one used to prove the IMS formula (19), we obtain (note that W{x) behaves like 
c|x| _1 and | VW(x)| behaves like c|x| -2 at infinity so there are no contributions 
from infinity to the integration by parts) 

j\VW 2 \ 2 = J \VF\ 2 \W\ 2 - J \F\ 2 WAW = j \VF\ 2 \W\ 2 

< sup \W(x)\ 2 [ r~ 2 dx = 4vrr sup |W(x)| 2 , 

|x|>r J2r<|x|<3r \ x \> r 

where the second equality follows since W is harmonic on the support of F. 
The estimate (105) now also follows from (104). 

We are now ready to estimate the TF densities. We shall use XtXrP tf 
for some R> r (possibly R is infinity) as a trial density in £ r 
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Since p° TF minimizes £° TF (p) + /x° TF / P and ^° TF = unless / p° T 
J xtp HF we have 

(107) /i° TF (| X r + /9 HP - | xtxRP TF ^j < £r TF (xtxRP T¥ ) ~ £r TF (P™] 

We write the right side as 
(108) 

sr (xtxnP^) - sr (pd = (xtx R p™) - sr (x+p tf ) 

+sr (x,V F )-^ OTF (P r ° TF )- 

We have for the first two terms 

(109) £ r ° TF (XrXRP™) - Sr" (xtp TF ) 

= J^ F P TF x + R + \\x + R P TF \\ 2 c 

< J V TF P TF X R + \\X R P TF \\ 2 c + crr-^' J X + R P TF , 
where we have used (104). For the last two terms in (108) we find 
(110) 

£r TF (XrP TF ) ~ £r TF (pr F ) = J W (p^ ~ Xr P™) ~ WP^ ~ Xr P 



,TF||2 

C 



+ / ( 

J\y\>r \ 



^ 2 ) 2/3 P TF (yf 3 -v TF (y)p TF (y) 



=^2)2/3^^5/3 _^ F(y)p OT F(y) 



dy. 



Using that p TF satisfies the TF equation (6) and that p TF = we see that for 
fixed y the expression 

^f'h^-^ F {y% t>0 

takes its minimal value for t = p TF (y) . Hence we conclude that the last integral 
above is negative. Thus combining (107)-(110) we have 



,OTF 



Mr" {J XrP HF ~ j XrXRP TF ^j < j W (p^ ~ xt P™) ~ {{P™ ~ xtp^Wl 

+ J ^ f p tf x + r + WxIp^WI + ot-^' J x + r p tf . 

Lemma 12.1 implies that J<p TF p TF X R , \\X R P TF \\ 2 C < (Const.)/?- 4 / X R P TF ■ We 
thus arrive at 
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(111) 

tf™ (/ *rP HF ~ j xtXnP 1 ^ < J W{P? T¥ ~ xtp TF ) - ||P r ° TF - XrVlc 

+ ((Const.)ir 4 + ar~ 4+£ ') J p TF X + R . 

This is the main estimate from which we shall derive the estimates of the 
lemma. We shall do this by choosing different values for R. One for the 
estimate (101) another for (102). We shall choose R such that 



(112) / Xr XRP TF < J XrP^. 

Let i? max denote the largest possible R for which this holds. Then 

(113) / X + R _P TF = (/ X,V F - / X,V F ) + < ar-^\ 

where the last inequality follows from Lemma 12.2. By Lemma 12.1 we have 
for all R > r > (3 Z~ 1/3 that 

(114) J XrP tf > (Const.)i?" 3 . 

We shall now make the assumption that D is chosen so small that if r < D 
then ar £ < 1. Thus from (113) and (114) we conclude that 

(115) R m 4 ax < (Const.)CT 4/3 r- 4+ t e ' < (Const.)r- 4 . 



From (113) and (111) with R = i? max we get, again using the above assumption 
on D, that 

(116) ||/f TF - X+P TF || C < J W{pP™ - XrP T¥ ) + (Const>r- 7+£ '. 
We estimate the integral on the right by dividing it in two parts 

J W(P? TF - xtp TF ) < j |^i|(p r OTF + XrP TF ) + j W 2 (P? TF - xtp TF ) 

< ar- A+e ' [ p° TF (x) + xt(x) P ' TF {x)dx + [ W 2 (p? TF - XrP TF ) 

Jr<\x\<3r J 

where we have also used (105). From Lemma 12.1 we arrive at 

(117) j W(p° TF - xU TF ) < (Const>r- 7 + £ ' + J VF 2 (p r OTF - xtp TF ). 

The last term in this estimate we now control using the Coulomb norm estimate 
Lemma 9.2. Note that p TF and p° TF both belong to L 6//5 since they are in 
L 5 / 3 n L 1 . We find from (106) that 

(118) J W 2 (p? TF - xtp TF ) < (Con S t.)c7r-i +£ '||p r OTF - Xr + P TF ||c- 
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Inserting the last two estimates into (116) gives (recall that ar £ < 1) 
||p r OTF " xtp TF \\ 2 c < (Const.)ar- 7+£ ' + (Const.) (ar^') 1/2 || p oxF _ x + p ^\\ Q 

and (101) follows immediately from this. 

We now return to (111) and make a new choice for R. It follows from the 
first inequality in (115) that that we can find R with R < R max satisfying 

iT 4 = (Const> 1 / 2 r- 4+£ '/2. 

Note that, since 

R- 4 < (Const.)a 1 / 2 ^'/2 r -4 ) 
we can choose D small enough to ensure that r < R. From (114) we have 

(119) J XrP tf > (Const.)a 3 / 8 r- 3+ f £ '. 

It follows from (113) and (114) that we may assume that the constant in the 
definition of R is chosen such as to ensure that / X^P TF > 2 J" X^ max P TF - Thus 
since (see (113)) 

/ *rV F " / XrXRP™ >-J X + R _P TF + / X + R P TF 

we have 

\ J X + R P TF < J XrP HF ~ J xtxRP TF . 

Thus from (119), (117), (118), (101), and (111) we conclude (102). □ 

Using these fairly weak estimates we shall now show that the outside TF 
potential and density satisfy Sommerfeld type estimates. 

Lemma 12.4 (Sommerfeld estimates for OTF). Assume N > Z. Given 
constants e',a > there exists a constant D > depending only on e',a such 
that for all r with ftoZ^ 1 ^ < r < D for which (85) holds for all \x\ < r then 
/i° TF = and for all \x\ > r we have 

(120) ^° TF 0c) < 3 4 2-V|xr 4 (l + ^r c |xr c ) 
and _ 2 

(121) tf TF (x) > 3 4 2- 3 7r 2 |x|- 4 (l + ar c , 

where a and A are universal constants (but not necessarily positive) and a > 
-1. Here as before ( = (-7 + v / 73)/2 « 0.77. 

Proof. Note first that the potential V r satisfies the assumptions in Theo- 
rem 4.6 with R = r. Hence if we can show that /x° TF < inf \ x \ =r ip° TF (x) the 
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potential y?° TF w hl satisfy the Sommerfeld estimates described in the theorem. 
In order to control tp° TF (x) for \x\ = r we note that for all \x\ > r we have 

(122) ^ r OTF (x) = <^ TF (x) + (x r V F - p r OTF ) * Ixl" 1 + W(x), 

where as in (103), W = <1 >HF — <3? TF . According to the Coulomb norm Corol- 
lary 9.3 we get for all s > that 



x^ tf -p° tf )*N _1 



< (Const.)s 1 / 5 max{||x^ TF || i 5/3 (B(;EiS)) ,||P r OTF |lL5/3(B( a; , s ))} 
+ (Const.) S - 1 /2|| Xr V F -p r OTF || c . 
From Lemma 12.1 we see that 

max{||XrV F || L 5/3 (B(a , ||P r ° TF |lL5/3(ij(., s ))} < (Const.)/- 6 * 9 / 5 , 

where we have assumed that D is such that ar 6 ' < 1. Inserting this and the 
estimate (101) from Lemma 12.3 above we obtain that for all \x\ > r 

(123) 

1/2 



XrV F -A° TF *M _1 



< (Const. )r~ 6 s 2 + (Const.)s _1/2 (crr" 7+£ ') 
= (Const> 2 / 5 r- 4+2£ '/5 



where we have optimized in s in order to get the last expression. Thus from 
(122), the iterative assumption (85), and Lemma 12.1 we obtain 

inf ^ r OTF (x) > inf v9 TF (x) - (Const.)a 2 / 5 r- 4+2e '/5 
\x\=r \x\=r 

> (Const.)/- 4 - (Const.)a 2 / 5 r- 4+2£ '/5. 

We have used that since ar e < 1 the error from (123) is worse than the error 
from (85). Note that the constant in front of r -4 above is positive. 

From Lemma 12.3 we know that n° TF < (Const. )cr 1 / 2 r~ 4+e '/ 2 Hence we 
may choose D such that n° TF < inf| x | =r ip° TF (x). 

As above we of course also have 



sup ip° TF (x) < (Const.)/- 4 + (Const.)a 2 / 5 r- 4+2e '/5. 

\x\=r 

Thus (Const.)/- 4 < inf ]x]=r <p° TF (x) < sup N=r ^ r OTF (x) < (Const.)/- 4 . 

,OTF 

V r we have 



That /U° TF = follows from Corollary 4.7 as follows. By harmonicity of 



lim (Ait)- 1 f r'V r {r'u)dw = (Air)- 1 f r$ HF (ro)dw = Z - [ XrP HF . 

r'->oo Js2 J S 2 r J 
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Thus from Corollary 4.7 we have if /x° TF / 0, i.e., if / p° TF = / Xr P HF that 
< ( / u r OTF ) 3 / 4 < (Const.) (Z-J XrP UF ~ J P r ° TF ) = (Const.)(Z - j p HF ), 

which is a contradiction since f p HF = N > Z. 

The estimates (120) and (121) now follow from Theorem 4.6. □ 

We are now ready to give the bounds on Ai and Ai defined in (98) 
and (99). 

Lemma 12.5 (Control of Ay and A2). Assume N > Z. Given constants 
e',a > there exists a constant D > depending only on e',a such that for 
all r with (3qZ^ 1 ^ < r < D for which (85) holds for all \x\ < r we have for all 
\x\ > r that 

(124) \Ai{r,x)\ < (Const.)r c |x|~ 4_c 
and 

(125) \A 2 {r,x)\ < (Const.)r c |x|- 4 - c , 
with C = (-7 + >/73)/2 w 0.77. 

Proof. Combining Theorems 5.2, 5.4, with Lemma 12.4 and recalling that 
fi TF = immediately gives the bound on A\. 

If we use the TF equation (33) we obtain from Lemma 12.4 that for all 
\y\ > r 

p? TF (y) - 3 5 2-V|y|- 6 | < (Const.)r^r 6 ^. 

Of course we similarly have from Theorems 5.2, 5.4, and the TF equation (6) 
that 

P TF (y) < 3 5 2'\\y\- 6 

and 

P TF (y) > 3 5 2- 3 vr|yr 6 - (Const.)Z-^ 3 |yr 6 - c . 
Since r > [5qZ~ 1 ^ we conclude that for all \y\ > r 

(126) \P TF (y)-P° TF (y)\ < (Const.)/|y|- 6 -^. 
Thus 

\A 2 1 < (Const.) / r C \y\~ 7 ^dy 

J\y\>\x\ 

which gives the bound in (125). □ 

We turn now to estimating ^,3. This requires estimating the difference 
between p° TF and p HF . We again begin by estimating this difference in the 
Coulomb norm. More precisely, we estimate in Coulomb norm the difference 
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between the "outside" TF density p° TF and the "outside" HF density Xr P™ '■ 
This is done through a semiclassical analysis of the exterior region {\x\ > r}. 

Lemma 12.6 (Coulomb norm comparison ofHFand OTF). Assume N > Z. 
Given constants e' , a > there exists a constant D > depending only on e',a 
such that for all r with (3qZ~ 1 ^ < r < D for which (85) holds for all \x\ < r 
we have 

(127) ||P° TF - xtp HF \\c < (Const.)r-i+^ 
and 

(128) J (x^P HF ) 5/3 < (Const.)r- 7 . 



Proof. Let 7 be the density matrix on L 2 (IR 3 ;C 2 ), but diagonal in spin, 
constructed in the semiclassical approximation Lemma 8.2 for the potential 
V = f° TF ■ Note that from Lemma 12.4 we have f° TF (y) > for \y\ > r and 
from its definition <p° TF (y) < for \y\ < r. From Lemma 8.2 we have that 



(129) 



p 7 = ^(stt 2 )- 1 [^r\T*9 



P° TF *9 2 , 



where we have used the TF equation (33) and the fact p° TF = proved in 
Lemma 12.4. Here g was given in Definition 8.1. From Lemma 8.2 we also 
have 



(130) 



Tr 



2 3/2 
5^" 



OTF] 5/2 



,1/2 r 



OTP]3/2 



(Note the factor of 2 in the formulas above compared to Lemma 8.2. This 
is due to the fact that there was no spin in Lemma 8.2.) The last inequality 
above follows from the Sommerfeld estimate for OTF given in Lemma 12.4. 

According to Lemma 8.5 we may, for all < A' < 1, choose a density 
matrix 7 such that its density p^ has support in {|x| > r} and such that 



P^< P~f and 



Tr 



(-±A-V r )7 



< Tr 



(-±A-V r )j\ + (Const.) / [V r ] 
+ (Const.)(AV)- 2 



x\<(;i-\')- 1 r 



M<(i 



5/2 
+ 



If we make the assumption that D is chosen to ensure that ar £ < 1 we may 
use the estimate on V r from Lemma 12.1 and use the same lemma to conclude 
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that p 7 (y) < (Const. )r 6 for all y. If we recall that V r has support for \x\ > r 
we get that 



Tr 



(-±A - V r ) 7] < Tr [(-3A - 7] + (Const.) (aV~ 7 + A'~ 2 r- 5 ) 



= Tr 



(-±A - y r ) 7] + (Const.)r~ 7+2/3 , 



where we have made the choice A' = r 2 / 3 and assumed that D is such that 
A' < 1/2. 

Since J P^ < J P 1 = J P r OTF < / xt P™ we see from Theorem 6.2 that in 
terms of the auxiliary functional £ A defined in (49) 

£ A [7 r HF ] < £A [7] + n < Tr [(-5A - F r ) 7] + (Const.)r~ 7+2 / 3 



+ 



3/7" P 7 (x)|x-y| 1 P 1 (y)dxdy + K, 



where 7 J HF and 7?. were defined in (48) and (51) respectively in terms of a 
parameter < A < 1 (different from the A' used above). We have here used 
that ^ F P7j = V r p^, since p^ has support in {\x\ > r}, and that 



Pz(x)\x-y\ 1 Pz(y)dxdy< \ \\ p 1 (x)\x - y\ 1 p 1 {y)dxdy 



1-1 



+ (Const. ) r - 7+2 / 3 



since p^ < p 1 . 

Since \x\~ x is super harmonic we have 

J J g(x — z) 2 \z — w\~ 1 g(y — w) 2 dz dw < \x — y\ 

and we conclude that 

£ A [7 r HF ] < Tr[(-lA-y r ) 7 

+ \fj Pr F (x)\x-y\- 1 pr F (y)dxdy + TZ 
From (129) and (130) we find that 

(131) £ A [7 r HF ] < ^° TF (P r OTF ) +Jv r (p r OTF - P r OTF * g 2 ) + (Const.) 
+ (Const.)r- 7+2 / 3 +K. 

We have 



s 2 r~ 3 



V r (p™-p™*g 2 



J (V r - V r * g 2 ) P ; 



2\ OTF 
r 



Now since V r (y) is harmonic for \y\ > r we conclude that V r * g 2 (y) = V r (y) 
for \y\ > r + s. Hence we get from Lemma 12.1 that 

(132) / V r ( Pr ° TF - /f TF * g 2 ) < (Const.)r- 4 / p r OTF < (Const.)r- 8 S . 

J J J\y\<r+s 
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We insert this into (131) and arrive at 
(133) 

£ A [7 r HF ] < £r° TF (P r OTF ) + (Const.) (s~ 2 r~ 3 + r~ 8 s) + (Const.)r~ 7+2 / 3 + K 
= £: r OTF (P r OTF ) + (Const.)r- 7+2 / 3 +^ 

with the choice s = r 5 / 3 . 

We shall now estimate 1Z. We shall choose the A used to define 7 ? HF and 
1Z in such a way that A < 1/2. We then see that the constant C\(r) in 
Theorem 6.2 satisfies C\(r) < (Const. )(Ar) -2 , since A < 1/2 and r < 1. We see 
from Lemma 10.5 that 

(134) / p HF < (Const.)r- 3 , 

J\y\>{\-X)r 

where we have used that r,ar £ ' < 1. 

Moreover, from Lemma 12.1 with r replaced by (1 — A)r we have 



/ [*(T-A)r (f )] T ^ ~ ( C ° nSt -) r ~ 7A - 



(l-A)r<| ?/ |<(l-A)-ir 

Hence we have 

11 < (Const.)A~ 2 r- 5 + (Const.)r^ 7 A + Ex [7 r HF ] ■ 

If we now use the exchange inequality in Theorem 6.4 and (134) we get (recall 
that pf F = # 2 p HF is the density corresponding to 7™) 

V2 , r _ 



^[7 r HF ] < (Const.) J (p^f 3 < (Const.) (|p r HF ) (/(O 
< (Const. )r~ 3/2 + r -7 ) 1/2 , 
where we have also used that according to (59) and Lemma 12.1 we have 

(135) J (p r HF ) 5/3 < (Const.)Te + (Const.)r- 7 . 
We may therefore conclude that 

(136) n < (Const.)r- 7 (r 2 A~ 2 + A) + (Const.)r^ 5 . 

We may use (135) and (136) to prove (128). Recall that P r HF (y) = P HF (y) if 
\y\ > (1 — A)~V. Now (128) follows if we simply observe that (135) and (136) 
hold with r replaced by r/2 and A = 1/2. We shall make a possible different 
choice of A below. 
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We shall now prove a lower bound on £ A [j^ F ] ■ We write 

£ A [7 r HF ] = Tr [(-±A - < F ) 7r HF ] +\JJ Pr(*)\x ~ Vl' 1 Pf" '(v)dx dy 
= Tr 



IA-C F + P° TF 



* \x\ 



1\ „,HP 



+ \\C-P° TF \\c 



p° TF (x)\x-y\- 1 p^(y)dxdy. 
If we use that on the support of p? F we have <&" F = V r we may write this as 



(137) 



£ A [7 r HF ] 



Tr 



+ ii/r-/T F ii 2 c 



(-^A-^ r OTF )7, 

pr F (^-y\- 1 Pr TF (y)dxdy. 



The trace may be bounded below by the sum of the first N' negative eigenvalues 
of the operator — ^A — y? 7 OTF , where N' is the smallest integer larger than 
Tr^™] = / pf F . From Lemma 8.2 (again with an extra factor of 2 due to 
spin) we therefore have that for all s > and all < 5 < 1 



(138) Tr[(-IA-^° TF )X P 



> 



2 5/2 



(157T 2 ) 

- 2Li<T 3/2 



l-5)- 3/2 /fe° TF ) 
P™ + 1 



5/2 



„OTP ,„OTF 
P r ~ % 



5/2 
5/2 



We first estimate the last term. Since p° TF * \x\ 1 is superharmonic we have 
by the mean value property that p° TF * |a;| _1 > P° TF * | rrr | 1 * g 2 . Thus we have 

9? r OTF - tf TF * g 2 = V r - V r * g 2 + p r OTF * l^r 1 * g 2 - p° TF * \x\- Y <V r -V r * g 2 . 

The same argument which led to (132) gives that V r (y) — V r * g 2 (y) = unless 
r — s < \y\ <r + s. Since by Lemma 12.1 we have |K(y)| < (Const. )r~ 4 (recall 
that V r is supported on {\y\ > r}) we obtain 



,„OTF ,„OTF , 2 



+ 



5/2 
5/2 



< (Const. )r~ s s, 



if we assume that s < r. (Note that s here does not have to be chosen as in 
the upper bound). 

From (134) we also get that / pf F < J p HF xt < (Const. )r~ 3 . 

Finally from the Sommerfeld estimate (120) and the fact tl 
positive only if V r (x) > 0, i.e., only if \x\ > r, we find that 

(<P° TF )T < (Const.)r- 7 . 



P? TF ( 



I- 
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We therefore see from (138) and the TF equation (33) (recall that ^° TF = 0, 
by Lemma 12.4) that if < 5 < 1/2 then 



Tr 



> ^(3vr 2 ) 2 / 3 J (p r OTP ) 5/3 - f ^ TF P° T 

-(Const.) (Sr~ 7 + <r 3 / 2 r _8 s) - (Const.)^ 2 ^ 3 
= U^ 2 f" J (^° TF ) 5/3 " / VT F /T F " (Const.)r- 7 +^ 

where we have chosen 5 = r ^ 2 / 5 s 2 / 5 an d s = r 11 / 6 , which agrees with s < r. 
If we insert this last estimate into (137) we obtain 

A r^/HFT ^ pOTF , QTF\ , II „HF _ „OTF||2 _ fr , n . \^~7+^ 



£ A [7 r HF ] > & (Pr TF ) + \\(>F - Pr TF \\c ~ (Const.)r 
If we compare this with (133) we see that 

\\Pr F ~ Pr TF \\c < (Const.)r- 7+ ^ + K. 

Finally, we use the Hardy-Littlewood-Sobolev inequality (80) and (128) to 
conclude that 



\Xr P — P r ||C 



< (Const.)||Xr + p HF - /H| 6 /5 < (Const.) ( [ P UF (yf /5 dy] 

/ \ 3/5 / \ 7/30 

< (Const.) / P^iyf'Hy) If ldy) 

\Jr<\y\ J \Jr<\y\<(l~\)-i-r / 



< (Const.)A 7 / 30 r- 7 / 2 . 
We thus get from (136) that 

IK/ + /iHF ^OTFII ^ II, +..HF ^HFII , ||„HF ^OTFII 

\\Xr P — P r \\CS\\XrP — P r Hc+llPj. - P r II G 
< (Const.)r-I+l + (Const.)r- 7 / 2 (rA- 1 + A 1/2 + A 7/30 ) 

which gives (127) if we choose A = min{l/2, r 5 / 7 }. □ 
We may now estimate A3 defined in (100). 

Lemma 12.7 (Controlling .A3). Assume N > Z. Given constants e',a>0 
there exists a constant D > depending only on e',a such that for all r with 
(3qZ~ 1 / 3 < r < D for which (85) holds for all \x\ < r we have for all \x\ > r 
that 

(139) \A 3 (r,x)\< (Const.)(|x|/r) 1/12 r- 4+ ^. 
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Proof. We shall use the Coulomb norm estimate (83) with / = /9? TF — 
X^P HF - We then immediately see from Lemma 12.6, and the fact, which follows 
from Lemma 12.4 and the TF equation (33), that p° TF (y) < (Const.)|y|- 6 , that 

\A 3 (r,x)\ < (Const.)(K|x|) 1/5 r~ 21/5 + (Const.) K~ l \x\~ l /2 r~ l z + l . 

This gives (139) if we choose k = (Const. )(r/|x|) 7 / 12 r 36 . □ 

End of proof of Lemma 10.4. For \x\ > f3$Z~ 1 / 3 the estimate in (87) 
follows from (128). For |x| < (3 Z~ 1/3 we get from (90) that 

/ P™(y) 5/3 dy < f p* ¥ {yf/ 3 dy < (Const.)Z 7 / 3 < (Const.)|x|~ 7 . □ 
J\y\>\x\ Jr* 

End of proof of the iterative Lemma 10.3. Let D > depending on a, e' 
be the smaller of the values D occurring in Lemmas 12.5 and 12.7. We may 
without loss of generality assume that D < 1. 

Given 5 > 0. We consider R < D satisfying /? Z~ 1/3 < Rl +S and we 
assume that (85) holds for all \x\ < Rq. 

Set R'o = Rl~ 5 and r = Rl +5 . Then we have (3 Z~ 1/3 < r < R < D and 
we can therefore apply Lemmas 12.5 and 12.7. Moreover R' Q > Rq. In order 
to prove (86) for Rq < \x\ < R' we use (97) and Lemmas 12.5 and 12.7. We 
obtain that for all \x\ > r 

4-1- J_ 



f 1 36 



$™(x) - *w(x)\ < (Const.) (r/|x|) C \x\~ 4 + (Const.) (|ar|/r) 
Moreover, for all Rq < \x\ < R' we have 

\x\ 1 - s < -r-r < \x\ d 

\x\ 

and thus 

$™\(x) - &w(x)\ < (Const.) (\x\- 4+<:5 + |x|" 4+ ^"^)^ . 

It follows that if 5 is small enough there exist e, C'$ > such that (86) is 
satisfied. □ 



13. Proving the main results Theorems 1.4, 1.5, 3.6, and 3.8 

The main result Theorem 3.6 on the maximal number of electrons iV is 
a simple consequence of Lemma 10.5 and Theorem 10.1. 

Proof of Theorem 3.6. We may of course assume that N > Z and that 
Z > 1 (otherwise the result follows from Lieb's bound Theorem 3.5). Then 
J p TF = Z. We can then use Lemma 10.5 with R chosen so small that Cm < 
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C$i? 4+e , because then (85) holds with a = 2C$ and e' = e. We conclude 
from Lemma 10.5 that 

[ p HF < f p TF (x) dx+aR- 3+£ ' + (Const.)(l+<7 J R e ')(l+ir 3 ) < Z+(Const.), 

J J\x\<R 

since now R, a, and e' are universal constants. We have thus concluded the 
result of Theorem 3.6. □ 

The asymptotics of the radius of an infinite atom given in Theorem 1.5 is 
a simple consequence of the main estimate Theorem 10.1 and the Sommerfeld 
asymptotics. 

Proof of Theorem 1.5. Note that in the neutral case N = Z we have from 
the main estimate Theorem 10.1 that 



I 



x\>R 



p TF (x)-p HF (x)dx 



[ P iiF {x)-p TF (x)dx 

J\x\<R 

(4tt)~ l R f ® T x(Ru>) - ^(Rco)duj 
J s 



< C^R- z+£ + C M R- 

Theorem 1.5 now easily follows from TF equation (6) and the Sommerfeld laws 
Theorem 5.2 and 5.4 for the case N = Z, i.e., p TF = 0. □ 

The potential estimate in Theorem 1.4 is somewhat more difficult to prove. 

Proof of Theorem 1.4. As in the proof of the main estimate Theorem 10.1 
we separately treat small \x\, intermediate \x\, and large \x\. 
We first consider small \x\. Note that 

^ F {x) - ^ F {x) = J (p TF (y) - p» F {y)) \x - y\- x dy. 

Thus using the Coulomb norm estimate (82) we obtain 

/ (P TF (y)-P HF (y))\x-y\- 1 dy 

< (Const.)s 1/5 max{||p TF || L 5/3 (B(X)a))) \\P UF \\ L ^(B(x,s))} 
+ (Const.) S - 1 / 2 ||(p TP -Ollc. 
If we use Lemma 11.1 and optimize in s we arrive at 

(140) \v UF {x) - 9? TF (x)| < (Const.)Z 4 / 3 - £3 < (Const.)|xr 4+3£3 - 4<5+3fea 

for some universal £3, if < I3qZ~ 1 ^. 

We now turn to intermediate \x\. We shall choose a D > such that, with 
the notation of Theorem 10.1, we have C M < C$D- 4+£ . Then by Theorem 10.1 
we have that (85) holds for all \x\ < D with a = 2C<j> and e' = e. We may now 
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assume that D < 1 and that D is smaller than the values for D in Lemmas 12.5 
and 12.6 corresponding to the above choices of a and e' . 
Consider 

((3 Z- 1/3 )TTs < | x | < I)TTs _ 

Set r = \x\ 1+s . Then \x\ > r and /?oZ -1 / 3 < r < D. We shall use the notation 
from Section 12. We write 



<p« F (x) - ^ F (x) = ^ F (x) - ^ OT » + ^° T » " 

The difference between the last two terms was defined in (98) to be Ai(r, x) 
and this was estimated in Lemma 12.5. 
We have 

^ F (x) - ^ TP (*) = J ( P ™(y) - xt(y)P HF (y)) \x - y^dy. 

Exactly as above, for small \x\, we now use Theorem 10.1, Lemma 12.6 and 
the Coulomb norm estimate (82) to conclude that 

\<p HF (x)-tf TF (x)\ < (Const.)r- 4+ 5T < (Const.)|x| (- 4 +^H 1+<5 ) , 

If we combine this with (124) from Lemma 12.5 we obtain 

(141) \ip nF (x) - ip TF (x)\ < (Const.) (|x|- 4+ ^ + |x|(- 4+ ^)( 1+5 )) . 

Combining (140) and (141) we see that by choosing 5 small enough we have 

i 

proved (10) for all \x\ < D 1 + s . 

i 

We turn to \x\ > D^ 6 , i.e., \x\ greater than some universal constant. 
Here we may write 



\^ F (x)-^ F (x)\< C(*)-C|C 



\ x \ K ' W 



x 



+ 



[ ( P TF (y)-P llF (y))\x-y\- 1 dy 

J\y\>\x\ 



The first term is controlled by the main estimate Theorem 10.1. If we use that 
according to Lemma 10.4 we have Ju^ui P™ {y) 5 ^ 3 dy < (Const.) and that the 
same estimate holds for the TF density (see Lemma 12.1) we may estimate the 
second term above as follows. Using Holder's inequality we have 



/ (P TF (y)-P HF (y))\x-yr 1 dy < (Const.) / \x - y^dy) 
•C , \J\x- y \<i / 

!- / (p TF (y)+P HF (y))dy< (Const.), 
J\x\<\y\ 



where the last estimate follows from Lemma 10.5. 



We end the paper by giving the proof of the bound on the ionization 
energy in Theorem 3.8. 
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Proof of Theorem 3.8. Since the HF energy is a nonincreasing function of 
N we have that < E HF (Z - 1, Z) - E UF (Z, Z). In order to prove an upper 
bound we shall construct a trial density matrix 7 for £ HF with Trpy] < Z — 1. 
We then clearly have that £ UF [7] > E UF (Z — 1,Z). Let be given in terms 
of appropriate r, A > as in the beginning of the proof of Theorem 6.2. We 
then choose as our trial matrix 

7 HF =e i H F0 _^ 

where 7 HF is the HF minimizer with Tr [7 HF ] = Z. According to the definition 
of we have 

Tr[ 7 « F ] <Z - { P HF (y)dy. 
J\y\>r 

We choose A = 1/2. Let R > be such that Cm = C<$>R~ i+e . We shall now 
choose r satisfying r < R. Then according to Theorem 10.1 we have that (85) 
holds for \x\ < r with a = 2C$ and e' = e. From Lemma 10.5 we therefore 
conclude that 



/ P™{y)dy = [p hf - [ P nF (y)dy 

J\y\>r J J\y\<r 

= I P TF (y)-p HF (y)dy+ f P TF (y)dy 

J\y\<r J\y\>r 

> [ P TF (v)dy - (Const. )r~ 
J\v\>r 



\<r J \y\>r 

-3+e 

l\y\>r 

where we have used that / p HF = J p TF = Z. 

We may of course assume that Z is larger than some fixed universal con- 
stant. For Z less than a universal constant, the total energy E HF (Z,Z) and 
hence the ionization energy E FiF {Z — 1, Z) — E FiF {Z, Z) are bounded by univer- 
sal constants (see Theorem 3.2). We can therefore assume that P Z~ 1/3 < R 
and we shall choose (3qZ^ 1 ^ < r. It then follows from Theorem 5.4, the 
TF equation (6) (recall that we consider the case p TF = and N = Z) that 
P TF {y) > (Const.) \y\- 6 for all \y\ > r. Hence 



/ p HF (y)dy > (Const. )r 3 - (Const.)r" 

J\v\>r 



\y\>r 

Thus we may choose r to be a small enough universal number (assuming that 
Z is large enough to allow PqZ^ 1 ^ < r) to ensure that f\ y \ >r P HF (y)dy > 1 and 

hence that Tr h HF j < Z - 1. 

From the estimate (53) in the proof of Theorem 6.2 we have 

£ HF [7 HF ] < £ HF [7 HF ] - £ A [7 r HF ] + ^ = ^ HF (^ Z) ~ £ A [% HF ] + ^ 

where as before 7 J HF = ^ r 7 HF ^ r , 1Z is given in (51), and the functional £ A was 
defined in (49). 
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It remains to prove that 
(142) -£ A [7 r HF ] + n < (Const.). 

As in (136) we conclude that 1Z < (Const. )r -7 < (Const.), where we have 
used that r is a universal constant. 

In order to estimate £ A [7 r HF ] we note that, since 3> r (y) is harmonic for 
\y\ > r and tends to at infinity, we have that for all \y\ > r 

< F (y) < sup <D r HF (x) < |y|-V sup \*™(x)\ + M- 1 (c $ r" 3 + £ + rC M ) 

\x\=r \x\=r 

< (Const. )r- 3 |y| _1 , 

where we have used the main estimate Theorem 10.1 and the bound on <3? TF in 

r 

Lemma 5.5 with fi TF = 0. If we use that r is some universal constant we get 

< F (y) < (Const.)lyr 1 . 

Using the Lieb-Thirring inequality (21) we see from the definition (49) of the 
auxiliary functional £ A that 

£ A [7 r HF ] > Kj P r(yf /3 dy - (Const.) / P -^dy 

^ ^ I y\ 

+ P r(x)\x-y\~ l pf F (y)dxdy. 

Here again p™ = 9^p HF is the density corresponding to 7 J HF . It follows from 
standard atomic TF theory that 

inf jtfi J p(y) 5 ^dy - (Const.) J ^dy+\ J f p{x)\x - y\~ l p(y)dx dy 

is some universal constant. Hence £ A [7^ F ] > —(Const.) and we have proved 
(142). □ 
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